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APRESENTAQÀO 


Este livro foi desenvolvido para dar suporte ao estudo de urna disciplina 
bàsica de Càlculo Numèrico. 

O Càlculo Numèrico tem por objetivo o estudo de métodos numéricos para 
a resolugào de problemas em geral e, em particular, destacamos neste texto 
os problemas ambientais. Os métodos numéricos serào desenvolvidos visando 
eficiència e estabilidade durante a aplicagào em problemas gerais. 

Consideramos de fundamental importància o Càlculo Numèrico no que 
se refere à aplicagào das disciplinas bàsicas corno Càlculo Diferencial Integrai, 
Geometria Analitica, Àlgebra Linear, Introdugào à Computagào, entre outras, e 
a utilizagào dos conceitos e resultados adquiridos nessas disciplinas, bem corno a 
aplicagào na resolugào numèrica de problemas. 

Apresentamos métodos numéricos para solugào de sistemas de equagòes 
lineares, equagòes algébricas, interpolalo e ajuste de fungòes, integragào nu¬ 
mèrica e equagòes diferenciais. 

Detalhamos alguns aspectos teóricos e computacionais, além da utiliza¬ 
gào de softwares especializados e reconhecidos na literatura, com a certeza de 
que estaremos contribuindo com ferramentas importantes para a elaboragào 
de projetos e pesquisas que estarào presentes nas atividades dos profissionais de 
cièncias em geral. 

Por firn, agradecemos a Carla Taviani Lucke da Silva Ghidini e Silvia Maria 
Pereira Grandi dos Santos pela colaboragào na revisào deste livro. 
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UNIDADE 1 


Teoria dos Erros em Processos Numéricos 



1.1 Introdugào 


Nas diversas àreas cientificas, diante da resolugào numèrica de um pro¬ 
blema, deparamo-nos com a necessidade de tratar com solugòes aproximadas, 
desde a interpretagào do problema, a modelagem matemàtica correspondente, 
erros nos dados obtidos por meio de medidas experimentais e implementagào 
de dados no computador. 

Para melhor entendimento, resumimos alguns passos considerados impor- 
tantes na resolugào de um problema da seguinte forma: 

1. Apresentagào do problema reai, entendimento e levantamento dos dados; 

2. Formulagào matemàtica correspondente do problema-modelo matemàtico; 

3. Elaboragào de um plano de resolugào e a escolha de um mètodo nu¬ 
mèrico adequado; 

4. Implementagào computacional do mètodo, linguagem de programagào; 

5. Anàlise sobre a coerència dos resultados obtidos e o problema inicial 
proposto; 

6. Reformular o modelo matemàtico e/ou escolher novo mètodo numèrico 
de resolugào, caso o Passo 5 nào esteja satisfeito. 


A execugào dos passos anteriores è chamada de Validagào do Modelo 
e as aproximagòes consideradas nesses passos levam a alguns tipos de erros, 
conforme exemplo a seguir: 


Exemplo 1.1 



Considere o problema de transporte de lixo contaminado, conforme exibido 
na Figura 1 : 


Figura 1 Lixo contaminado. 






Urna transportadora possui très tipos de caminhòes representados por C r 
C 2 e C 3 , os quais sào equipados para levar très tipos diferentes de materiais con- 
taminados M i; M 2 e M 3 para lixòes adequados, conforme a Tabela 1.1: 


Tabela 1.1 



MI 

M2 

M3 

CI 

6.8 

1.9 

1.7 

C2 

1.5 

2.6 

1.4 

C3 

2.3 

1.8 

2.3 


O caminhào C 1 pode levar 6.8 toneladas do tipo M 15 1.9 tonelada do tipo M 2 
e 1.7 tonelada do M 3 . Deseja-se saber quantos caminhòes de cada tipo devemos 
enviar para transportar 16 toneladas do tipo M r 7 toneladas do M 2 e 6 toneladas 
do M 3 ? 

Exibiremos a execugào de todos os passos citados anteriormente e fare- 
mos urna anàlise final dos resultados obtidos. 

• Passo 1 

Neste passo, teremos o entendimento do problema observando todos os 
detalhes, caracteristicas e condigòes intrinsecas do mesmo. 

• Passo 2 

Faremos a formulagào matemàtica do problema, para isso, devemos deter¬ 
minar as variàveis do problema. 

Assim, seja o nùmero de caminhòes x jt j = 1,2,3 a ser determinado para 
o transporte do lixo contaminado. Dessa forma, observando a tabela de dados 
vemos que podemos escrever a seguinte equagào para o transporte do lixo con¬ 
taminado 

6.8 x 1 +1.5 x 2 +2.3 x 3 =16 

Seguindo o mesmo raciocinio, para o transporte dos lixos contaminados 
M 2 e M 3 , temos o seguinte sistema de equagòes lineares para ser resolvido, que 
sera estudado nas próximas unidades de estudo. 

6.8x 1 + 1.5 x 2 + 2.3x 3 = 16 
< 1.9x 1 + 2.6 x 2 + 1.8x 3 = 7 
1.7x 1 + 1.4 x 2 + 2.3x 3 = 6 
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Passo 3 


Devemos escolher o mètodo adequado de resolugào. Neste caso, vamos 
utilizar o Mètodo de Eliminagào de Gauss, o qual sera alvo de estudo nas pró- 
ximas unidades. 

• Passo 4 

Neste passo, podemos implementar no computador o mètodo de resolugào 
escolhido ou utilizar um software apropriado e teremos a solugào do sistema 
corno segue: 


x 1 = 1.9452, x 2 = 0.7953 e x 3 = 0.6868 


• Passo 5 

Temos que analisar a solugào encontrada, pois o nùmero de caminhòes 
deverà ser um nùmero inteiro positivo. Como temos urna solugào aproximada 
nào inteira para o sistema de equagòes lineares, podemos aproximà-la por um 
nùmero inteiro mais próximo de tal solugào. Dessa forma, tomamos corno solu¬ 
gào x 1 =2.0, x 2 =1.0 e x 3 = 1.0. 

Assim, devemos contratar dois caminhòes do tipo C i; um caminhào do tipo 
C 2 e um caminhào do tipo C 3 . 

Observagòes 

a) Na modelagem matemàtica do problema, consideramos urna aproxima- 
gào para o problema dado, urna vez que a capacidade de transporte 
de cada material em cada caminhào è estimada, isto è, com urna certa 
margem de erro. 

Alguns detalhes foram desconsiderados, corno a integralidade da solugào 
e custos envolvidos no processo de transporte. 

b) A solugào obtida durante a execugào do mètodo de resolugào foi apre- 
sentada com urna aproximagào de quatro casas decimais, o que signifi¬ 
ca que estamos cometendo erros em todas essas aproximagòes. 

Diante das consideragòes anteriores sobre o tratamento dos problemas com 
aproximagòes, desde a modelagem matemàtica, os erros cometidos durante o 
processamento dos métodos de resolugào no computador, è necessàrio um en- 
tendimento geral sobre tipos de erros existentes, corno descreveremos a seguir. 




1.2 Erros na fase da modelagem 


Devido às simplificagòes no processo de modelagem matemàtica de um 
problema, que muitas vezes sào necessàrias, podem ocorrer erros na represen- 
tagào do fenòmeno da natureza que estivermos analisando. 

Os problemas ambientais geralmente sào complexos e o modelo matemà¬ 
tico é urna aproximagào do problema reai, representado por expressòes mate- 
màticas que muitas vezes necessitam de algumas simplificagòes para obtermos 
urna solugào aproximada que resultam em erros. 

1.3 Erros na fase de resolugào 

Sào erros devido ao fato dos equipamentos computacionais terem capacida- 
de limitada para armazenar os dfgitos significativos de valores numéricos, utilizados 
nas operagòes elementares de adigào, multiplicagào, subtragào, divisào, etc. 

1.4 Erros na mudanga da base 

Considerando que a maioria dos computadores representa os valores nu¬ 
méricos na base binària quando sào armazenados, estes sào transformados 
em geral, da base decimai para a base binària ou em outra representagào. Essa 
transformagào pode ser acometida de erros, devido à limitagào da representa- 
gào do equipamento computacional que estamos utilizando para o processa- 
mento dos dados numéricos. 

Dado um nùmero reai N, é sempre possivel representà-lo em qualquer 
base b, da seguinte forma: 


N > : = Z a i b ' 


em que a, e {0,1,2,3,...,(b-1 )}, com n e m inteiros. 

Exemplo 1.2 

Numeros escritos na base 2 

a) (1011) 2 =1*2° + 1 *2 1 + 0*2 2 + 1*2 3 

b) (111.01) 2 = 1 *2 2 +0*2~ 1 + 1*2 0 + U2 1 + 1 *2 2 




Exemplo 1.3 

Numeros escritos na base decimai 

a) (231 ) 10 =1*10° + 3*10 1 + 2*10 2 

b) (231.35) 10 = 5*10 2 + 3*10' 1 + 1*10° + 3*10 1 + 2*10 2 

Mudanga da base decimai para a base binària 

Nùmero na base decimai com somente a parte inteira. 
Procedimento padrào: Divisòes sucessivas. 

Exemplo 1.4 

Podemos escrever o nùmero (25) 10 na base 2, corno segue: 

(25) 10 = 1*2° + 0*2 1 + 0*2 2 + 1*2 3 + 1*2 4 
= ( 11001) 2 

isto é, aplicando o processo de divisòes sucessivas temos: 

25+ 2 = 12 e resto = 1 
12+2=6 e resto = 1 
12+2 = 6 e resto = 0 
6+2 = 3 e resto = 0 
3 + 2 = 1 e resto = 1 


Podemos esquematizar o processo de divisòes sucessivas, usando o Mètodo 
da Chave, da seguinte forma: 


25 

1 


12 

0 


2 

6 

0 


2 

3 

1 


2 

1 

1 


2 

0 


Para obter o nùmero 25 na base 2, basta tornar os restos das divisòes su¬ 
cessivas, de 25 por 2, do quociente 12 por 2, e assim por diante, do final para o 
infoio no esquema anterior. 


Assim, escrevemos o nùmero 25 corno segue: 


25 = 1*2° + 0* 2 1 + 0*2 2 + 1*2 3 + 1*2 4 
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E, portanto temos: 


(25) 10 ={11001) 2 

Definicào 1.1 Erro de Arredondamento 

Dizemos que um nùmero x foi arredondado na posigào n se todos os digitos 
de ordem (n + 1) sào desprezados da seguinte forma: 

O digito de ordem n é acrescido de urna unidade se o de ordem (n + 1) for 
maior ou igual a 5. Caso contràrio, o nùmero x é representado apenas com os 
n digitos iniciais. 

Exemplo 1.5 

a) O nùmero x = 0.123577, usando a regra de arredondamento anterior, obte- 
mos o nùmero x arredondado com 4 casas decimais da seguinte forma: 


x' = 0.1236 


b)0 nùmero x = 123343, usando a regra de arredondamento anterior, obte- 
mos o nùmero x arredondado com 3 casas decimais da seguinte forma: 


x' = 0.123 


Definigào 1.2 Erro deTruncamento 

Quando representamos urna fungào através de urna sèrie infinita, o erro no 
valor de f{x) ao truncarmos a sèrie após um nùmero finito de termos è chamado 
de erro de truncamento. 

Exemplo 1.6 

a)Consideremos a representagào de urna fungào f(x) utilizando a Sèrie 
de Taylor, nas vizinhangas do ponto x, corno segue: 

m = m + + r g >(x) (x g f ) %... + f<" > (x) (x n f ) " +... 


em que, f (n \x) è o valor da n-ésima derivada da fungào f{x) no ponto x. 

Quando truncamos a sèrie no 3 a termo, isto è, considerando apenas os 
termos até a derivada de ordem 2, na expressào acima, temos um erro cometido 
nessa aproximagào, corno segue: 









b)Consideremosodesenvolvimentodafungào f{x) = e x em Sèrie de Taylor, 
nas vizinhangas do ponto x = 0, isto é: 

Como, f{x) = e x , f (1 \x) = e x , f (2 \x)e x ... 

Ainda, f(0) = 1, f (f \0) = 1, f (2) (0) = 1... temosque: 


e x = 1 + x + — + — + ... + — + ... 


2! 3! 


n! 


Considerando apenas os quatro primeiros termos da sèrie, temos: 


y2 y 3 i 

e x = 1 + x + — + — =- (x 3 + 3x 2 + 6x + 6) 
2! 3! 6 


Para x = 2, temos e 2 = 6.33333, que è um valor com erro absoluto bem 
significativo, quando comparado com o valor e 2 = 7.38906 obtido numa calcu- 
ladora cientffica que armazena urna quantidade maior de termos da sèrie. 

Definigào 1.3 Erro Absoluto 

Quando substituimos urna solugào exata x de um problema por urna solugào 
aproximada x', cometemos um erro chamado Erro Absoluto: 



Definigào 1.4 Erro Relativo 

Quando consideramos o Erro Absoluto cometido em relagào a urna gran- 
deza numèrica, chamamos de Erro Relativo: 








Exemplo 1.7 


a) Consideremos o valor exato x = 2345.713 e o valor aproximado 
x = 2345.000. 

Nesse caso temos: 

E a = | 2345.713 - 2345.000 \ = 0.713 e, 

E = = 0713 = 0.00030396 

r | x | 2345.713 

b) Consideremos o valor exato x = 1.713 e o valor aproximado x' = 1.000. 


Entào, 


E a =\1.713-1.000\ = 0.713 

E r = = °^ 3 - = 0.416229 

r | x | 1.713 

Observe que nos exemplos anteriores a) e b) o erro absoluto é o mesmo, 
embora o erro cometido pela aproximagào seja muito mais significativo no exem¬ 
plo b). No exemplo a), o erro relativo é da ordem de 0.03% enquanto no exemplo 
b) é da ordem de 41,6%. 

Exemplo 1.8 

Considere as distàncias entre as cidades: 

a) Sào Carlos a Campinas -170 km 

b) Sào Carlos a Ibaté - 20 km 


Ao ser perguntado sobre qual a distància no item a) e b), urna pessoa in¬ 
forma corno a seguir: 

a) Sào Carlos a Campinas -160 km 

b) Sào Carlos a Ibaté -10 km 








Qual o erro cometido nessas informagòes? 


E a =\l70 -160 \ = 10 ema) 


E a = \20 -10\ = 10 emb) 



10 

170 


em a) 



10 _ 

20 


1 

— emb 
2 


Podemos concluir que o Erro Absoluto em a) ou b) é o mesmo, nào for- 
necendo informagòes significativas sobre o quanto estamos errando naquela 
informagào. 

O Erro Relativo nos informa que ao errar 10 em 170 nào foi cometido um erro 
muito grande na informaqào, pois a grandeza numèrica 10/170 = 0.0588 = 5.88%, 
mas ao errar 10 em 20, foi cometido um erro de 10/20 = 0.5 = 50% na informagào. 

Podemos concluir que o Erro Relativo nos fornece mais informagòes sobre 
a qualidade do erro que estamos cometendo num determinado càlculo, urna vez 
que no Erro Absoluto nào é levada em consideragào a ordem de grandeza do 
valor calculado, enquanto que no Erro Relativo essa ordem é contemplada. 

1.5 Erro Absoluto e Erro Relativo nos procedimentos numéricos 

Como utilizar as definiqòes de Erro Absoluto e Erro Relativo? 

Sabemos que nos procedimentos numéricos em geral, quando resol- 
vemos um problema, geramos urna sequència de soluqòes aproximadas 

^0’ ,■■■ ,X k , Xk+1’■■■’ X- 

Se essa sequència for convergente para a solugào x, na medida em que k 
cresce, k °° as soluqòes aproximadas x k e x k+1 tendem a ficar próximas. 

As definigòes de E a e E r definidas anteriormente serào usadas para medir 
o quào próximas estào às solugòes x k e x k+1 e dessa forma podemos inter¬ 
romper a sequència x k gerada, utilizando um Critèrio de Parada, com o Erro 
Absoluto ou o Erro Relativo, corno segue: 







Dada urna tolerància por um nùmero £ > 0, suficientemente pequeno, con- 
sideramos que: 


E a = \ x k+1 -x k |< e 


E = ' Xk+1 * k I < £ 
x k+1 


Caso E a ou E r seja menor do que o nùmero e, entendemos que as solugòes 
x k e x k+1 estào próximas e podemos interromper a sequència tornando o ùltimo 
valor calculado corno a solugào aproximada do problema. 

Exemplo 1.9 


Podemos resolver a equagào x -2 = 0 utilizando o processo iterati¬ 
la 
2 


vo x k+1 = - 


( 2 ' 

x k + — 


V 


,k = 1,2,.... A medida que variamos k = 1,2,... gera- 


y 


mos urna sequència de solugòes aproximadas x 0 , x u ... que converge para 
x = 42 = 1.414213562. 

Para urna tolerància dada e = 0.0001 e x 0 =1, urna solugào inicial, temos 
a sequència de solugòes aproximadas convergente para a solugào x = 42, ge- 
rada pelo processo iterativo, isto è: 

Tornando x 0 = 1, urna solugào inicial qualquer e, no processo iterativo dado, 
variando k = 1, 2,... temos: 


Para k = 0^>x 1 = 1.50000 




= 0.33333 > e 


Para k = 1^>x 2 =1.41667 


x 2 ~x , 


= 0.05882 >e 


Para k = 2 ->x 3 = 1.41422 -> ^ * 2 ' =0.00173 >e 


Para k = 3->x 4 = 1.41421 -» ^ X . 3 ' =0.00001 < e 


'4 


22 











Dessa forma, usando o Erro Relativo corno o Critèrio de Parada, e urna 
vez que este està satisfeito na sequència quando k = 3, temos que a solugào 
x = x 4 =1.41421. 

Nota: 

Podemos ainda considerar os erros de representagào que, em geral, sào 
utilizados nos computadores para representagào dos numeros, chamados de 
ponto flutuante. 

A uniào de todos os numeros na representagào em ponto flutuante é cha- 
mado de Sistema de Ponto Flutuante. Detalhes sobre este tòpico podem ser 
vistos na Referència [1]. 

1.6 Propagagào dos erros 

A propagagào dos erros pode ser observada quando utilizamos um pro¬ 
cesso numèrico para buscar a solugào de um determinado problema. Esse 
processamento envolve um nùmero muito grande de operagòes elementares, e 
os erros acumulados no final dos processos numéricos interferem na qualidade 
da solugào aproximada do problema. 

1.7 Mapa Conceitual 

Um Mapa Conceitual è um instrumento utilizado no processo de ensino, 
aprendizagem e avaliagào, por meio de urna representagào gràfica de um con- 
junto de significados conceituais, o qual pode vir acompanhado de um pequeno 
texto explicativo e de simbolos que ajudam a organizar, sistematizar, estudar e 
detectar as ideias principais do tòpico abordado. 

Por exemplo, apresentamos a seguir um Mapa Conceitual preliminar sobre 
os erros em processos numéricos. 





I 

t 

f 1 

Sotucào Aproximada 

v_!_/ 

Figura 1.2 Um Mapa Conceitual sobre erros. 
































1.8 Exercfcios 


1. Converta os numeros decimais para a base binària: 
a) 45 b) 2978 c) 7699 


2.Resolver a equagào x 2 -7 = 0 usando o seguinte processo iterativo: 


X k+1 ~ ' 


1 


x k +- 


, k = 0,1,... com e =10 4 no critèrio de parada partindo de 


um ponto inicial x 0 = 1. 

3. Usando a tècnica de arredondamento, representar os numeros abaixo 
com 4 casas decimais: 


a) 0.4567897 b) 67.456789 c) 12.677777 


4.Sabemos hoje que a distància mèdia da Terra à Lua è de 384000 km. A 
distància da Terra à Lua obtida por cientistas da antiguidade, na segunda meta- 
de do séc. Il a.C., foi de 402500 km. 


a) Calcule o erro absoluto cometido nessas informagòes. 

b) Calcule o erro relativo cometido. Qual è o valor do erro relativo em 
porcentagem? 

5.Faga um mapa conceitual detalhado sobre os erros em processos nu- 
méricos, considerando as definigòes, introduzindo labels indicando se aprendeu 
(A), nào aprendeu (N), se gostou (G), detestou (D ), achou interessante (I), etc. 

Observagào: Os mapas conceituais podem ser feitos com os softwares dis- 
poniveis corno: Xmind <http://www.xmind.net>, CMAPS <http://cmap.ihmc.us> 
ou mesmo a barra de desenhos do Office, Word ou Powerpoint. 







UNIDADE 2 


Solugào Numèrica de Sistemas 
de Equagòes Lineares 



2.1 Introdugào 


O conhecimento de métodos numéricos para resolugào de sistemas de 
equagòes lineares torna-se necessàrio devido ao fato de que muitos problemas 
após sua modelagem matemàtica correspondente resultam na resolugào de um 
sistema de equagòes lineares. 

Esses sistemas lineares em geral sào de grande porte, a matriz A dos coe- 
ficientes possui muitas linhas e colunas, e necessitam de métodos numéricos 
estàveis, isto é, métodos em que os erros de arredondamento nào destroem a 
qualidade da solugào aproximada obtida, e com pouco tempo de processamen- 
to nos computadores. 

Os sistemas de equagòes lineares podem ser quadrados, quando o nù¬ 
mero de linhas é igual ao nùmero de colunas, ou retangulares, caso contràrio. 
Denominamos sistemas consistentes quando apresentam pelo menos urna so- 
lugào, e inconsistentes, caso contràrio. 

Para ilustrar os sistemas de equagòes lineares, considere o problema de 
composigào de areias para filtros de purificagào de àgua, conforme Figura 2.1 
a seguir: 



Figura 2.1 Filtros de purificagào de àgua. 

Devido a sua propriedade de grande permeabilidade, a areia é usada na cons- 
tituigào de filtros de Estagòes de Tratamento de Àguas de abastecimento (ETÀ) 
corno meio filtrante, por interceptar as impurezas existentes na àgua afluente. 

Para construgào das unidades de filtragào de urna ETÀ, urna usina de agre- 
gados foi contratada para fornecer a areia destinada à filtragào. Para atender ao 
pedido, a usina dispòe de areias brutas provenientes de très portos passiveis de 
exploragào com composigòes granulométricas distintas, conforme Tabela 2.1 : 









































Tabela 2.1 


Faixa granulométrica (mm) 

Porto 1 

Porto 2 

Porto 3 

0.42-0.59 

0.28 

0.10 

0.15 

0.59-0.71 

0.12 

0.22 

0.10 

0.71-0.84 

0.14 

0.12 

0.18 


Na Tabela 2.1, tornando 1m 3 da areia bruta proveniente do Porto 1 tèm-se 
0.28m 3 de gràos na faixa 0.42-0.59mm. As areias brutas sào peneiradas para 
obter as areias nas faixas especificadas. 

Essas areias serào dispostas em camadas que deverào obedecer às com- 
posigòes granulométricas estabelecidas por norma, nas quantidades a seguir, 
conforme Tabela 2.2: 


Tabela 2.2 


Faixa granulométrica (mm) 

Volume de areia (m 3 ) 

0.42-0.59 

20 

0.59-0.71 

10 

0.71-0.84 

12 


Deseja-se determinar quais as quantidades de areia de cada porto devem 
ser extraidas de forma que atendam a demanda. 

Na modelagem matemàtica do problema, sejam x y , j = 1,2,3 as quantida¬ 
des de areia em m 3 a serem retiradas dos respectivos portos P t , i = 1,2,3. 

Dessa forma, para obtermos as quantidades de areia de cada porto que 
devem ser extraidas de forma que atendam a demanda, temos que encontrar a 
solugào do seguirne sistema de equagòes lineares: 

0.28x +0.10x +0.15x =20 

1 2 3 

< 0.12x +0.22x + 0.1 Ox =10 

1 2 3 

0.14x +0.12x +0.18x =12 

li 2 3 

Observagòes 

a) Nesse problema algumas simplificagòes foram feitas, por exemplo, nào 
foram considerados os custos de extragào de areia nos portos, e ainda 
consideramos apenas très portos de extragào. À medida que aumentamos 
o nùmero de portos, os sistemas terào mais equagòes e mais variàveis. 

b) A resolupào deste sistema de equagòes, requer o conhecimento de mé- 
todos adequados e eficientes de resolugào, quais serào apresentados 
nesta unidade, corno segue: 

















Considere o sistema de equagòes lineares Ax = b, em que A=(a jj ), 
i,j = 1,...,n é matriz dos coeficientes, x = (x p x 2 ,xj o vetor incògnita 
transposto, e b = (b p b 2 ,...,b n Y é o vetor dos termos independentes transposto. 
Consideramos que o determinante de A, det( A) * 0, o que nos garante a unicidade 
da solugào para o sistema linear. 

Representamos o sistema de equagòes lineares da seguirne forma: 


a x + a x + ... + a x = b 

11 1 12 2 In n 1 

a x + a x + ... + a x = b 

12 1 22 2 2n n 2 


a x + a x +... + a x 

ni 1 n2 2 nn n 


b 

n 


Na forma matricial temos: 


(a a ... a 'ì 


fx ì 


fb } 

11 12 In 


1 


i 

a a ... a 


X 


b 

21 22 2n 


2 

= 

2 

a a ... a 


X 


b 

^ ni n2 nn ^ 


K n ) 


k n y 


A x b 


Resolver o sistema anterior dado consiste em determinar um vetor 
x = (x u x 2 ,...,x n ) t que satisfaga todas as equagòes simultaneamente. 

Podemos representar graficamente a solugào de um sistema no 3ì 2 , corno 
no exemplo a seguir: 

Exemplo 2.1 

x + x =2 

1 2 

< 

2x + x =3 
l 1 2 

Observe que det( A) = -1*0 e que a solugào x = (1,1) 'encontra-se na in- 
tersecgào das duas retas, representando as duas equagòes lineares, conforme 
Figura 2.2: 












Graficamente temos: 



Figura 2.2 Solugào de um sistema linear no plano. 

Podemos observar na Figura 2.2 que a solugào do sistema de equagòes 
encontra-se na intersecgào das retas, isto é, no ponto x 

Podemos ainda representar a solugào de um sistema de equagòes no es- 
pago Sì 3 , isto é, um sistema de equagòes com 3 equagòes e 3 incógnitas, con¬ 
forme a seguir: 

Exemplo 2.2 


2x 

+ 

X 

+ Ox = 

3 

i 


2 

3 


-4x 

+ 

2x 

+ X = 

-2 

i 


2 

3 


Ox 

i 

+ 

4x 

2 

+ X = 

3 

4 


Nesse caso, a solugào x = (1 , 1,0)' do sistema de equagòes encontra-se 
na intersecgào dos très planos representados pelas equagòes do sistema dado, 
conforme ilustrado na Figura 2.3 a seguir: 
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Figura 2.3 Solugào de um sistema linear no espago. 








Métodos de Resolugào 

Na resolugào de sistemas de equagoes lineares, podemos utilizar métodos 
diretos ou métodos iterativos. 

2.2 Métodos diretos 

Um mètodo direto para resolver o sistema Ax = b, em que A= (3//). 
i, j = 1n com det(A)*0, consiste em determinar exatamente o vetor solugào 
x = (x 1 ,x 2 ,...,x n ) t , se nào fossem os erros provenientes do processamento de 
càlculos envolvidos nos algoritmos. 

2.2.1 Sistema triangular inferior 

Seja Ax = b um sistema de equagòes lineares, em que A= (*//). i,j = 1,—,n 
é triangular inferior, isto é, os seus coeficientes (a /y ) = 0 sempre que i < j e com 
a n *0 para 




Representamos um sistema triangular inferior por: 


a x 

ii i 


a x + a x 

21 1 22 2 



a x + a x +... + a x = b 

_ ni 1 n2 2 nn n n 


Para construir o algoritmo para resolugào do sistema triangular inferior, 
destacamos a linha genèrica (/'), isto é: 

a n X 1 + a, 2 x 2 +... + a (/ x ( - =bj 


Considerando a /( * 0 e isolando a incògnita x ( na equagào anterior, temos: 


x ( . = (b-, - a i1 x 1 + a i2 x 2 + a i3 x 3 +... + a ff _,x,._,) /a n 





ou seja, 


^ij x i 
i=i 


Dessa forma, podemos escrever o algoritmo para resolugào de sistema 
triangular interior corno segue: 

Algoritmo 2.1 

Faga: 



Para / = 2,...,n, faga: 

i-1 

1 a ii x i 

j = i 


Exemplo 2.3 

Considere o sistema de equagòes lineares triangular interior: 


~6 

0 

0~ 

V 


~6 

1 

5 

0 

x 2 

= 

6 

0 

1 

1 

_ X 3_ 


1 


Usando o Algoritmo 2.1 temos: 



_ b 2 -a 21 x 1 _ 6-1(1) _ 1 
3^2 5 

_ b 3 - a 31 x, - a 32 x 2 _ 1 - 0(1) - 1(1) _ Q 

a 33 



1 
















Portanto, temos a solugào do sistema x = ( 1,1, 0) f . 


2.2.2 Sistema triangular superior 


Seja Ax = b um sistema de equagòes lineares, em que A={a jj ), /', j = 1,...,n 
é triangular superior, isto é, os seus coeficientes (a, 7 ) = 0 sempre que i> j e 
com a /; * 0, para 

i = 


Representamos um sistema de equagòes triangular superior por: 
a„x, + a 12 x 2 + ... + a 1n x n = b, 

a 22 X 2 + ■ ■ ■ + a 2n X n = 2 

< 

a nn X n ~ ^n 


Para a construgào do algoritmo destacamos a linha genèrica (/), isto é: 


a ii X i + a a + l X i+1 + a a + 2 X i+2 + a ii+3 X i+3 + ■■■ + a in X n - b i 


Considerando a„ * 0 e isolando a incògnita x, na equagào anterior, temos: 


n 

I a n x i 

j=i+1 


Dessa forma, podemos escrever o algoritmo para resolugào de sistema 
triangular superior da seguinte forma: 

Algoritmo 2.2 

Faga: 



x 


n 







Para / = faga: 
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X a u x j 

i= i+i 


Exemplo 2.4 

Considere o sistema de equagòes lineares triangular superiori 


5 

1 

0~ 

x i 


6 

0 

7 

-1 

X 2 

= 

6 

0 

0 

8 

_ X 3_ 


8 


Usando o Algoritmo 2.2 temos: 



_ b 2 &23 X 3 _ 6 ( 7)(7) _ 
2 _ ~ -, 



x _ b, - - a 73 x 3 _ 6 - (7)(7) - (0)(7) _ ? 

? a„ 5 

Portanto, temos a solugào do sistema x = {1,1,1)'. 

Observagào 

O Esforgo Computacional E c de um Algoritmo é a quantidade de opera- 
gòes elementares envolvidas durante a aplicagào do mesmo. 

No caso da solugào de um sistema de equagòes triangular superior ou interior, 
o esforgo computacional dos Algoritmos 2.1 e 2.2 é E c = rt 2 operagòes elementa¬ 
res, sendo ri o nùmero de operagòes de divisào, n ^ n ^ o nùmero de operagòes 
de adigào (ou subtragào) e n ^ n — o nùmero de operagòes de multiplicagào. 

No exemplo anterior, temos queo esforgo computacional E c = n 2 ={3) 2 =9, 
em que très operagòes sào de divisào, très operagòes sào de adigào ou subtra¬ 
gào e très operagòes sào de multiplicagào. 

















Experièncias computacionais mostram que o tempo computacional envol- 
vido nessas operagòes é pequeno, tornando os sistemas triangulares bastante 
atrativos, os quais serào usados no mètodo de decomposigào posteriormente. 

Definigào 2.1 Menores principais 

Seja a matriz A = (a ij ), i,j = 1,...,n. Denominamos menores principais de 
ordem k da matriz A por: 

A k =det{A k ) 


em que A k =(a /y ), i,j = 1,...,k è formada pelas k primeiras linhas e k primeiras 
colunas de matriz A. 

Exemplo 2.5 

Considere a matriz: 


A = 


5 

1 

1 


0 

2 

2 


0 

0 

1 


Càlculo dos menores principais: 


A, 






0 

2 

2 


0 

0 

1 


= 10 


Métodos de Decomposigào 

Os métodos de decomposigào para resolugào dos sistemas de equagòes 
lineares Ax = b consistem em decompor de forma ùnica a matriz A do sistema 
no produto de duas matrizes triangular interior e superior respectivamente, e 
resolvermos dois sistemas triangulares interior e superior respectivamente, os 
quais sào facilmente resolvidos corno os exibidos anteriormente. 


2.2.3 Mètodo de Decomposigào LU 


O Mètodo de Decomposigào LU consiste na decomposigào ùnica da matriz 
A = (a jj ), i,j = do sistema de equagòes, no produto de urna matriz L = (£%), 
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i,j = 1,...,n triangular interior, com os elementos da diagonal £ ti =7, / = 1,..., n por 
urna matriz U = ( u /y -), i,j = 1,n triangular superior. 

Dessa forma, enunciamos 0 seguirne resultado: 

Teorema 2.1 

Seja urna matriz A=(a ij ), i,j = Seos menores principais de A, A, *0, 
i = 1,2, n — 1, entào A se decompòe, de maneira ùnica, no produto de urna ma¬ 
triz triangular interior L=(f ij ), i,j = com £ r ,=1, por urna matriz triangular 

superior U={u ij ),i,j = 

n 

Além disso, det (A) = det (U) = Y[ u a 

i = 1 

Prova: Referència [1]. 

Processo de Decomposigào LU 

Por facilidade de entendimento, vamos decompor urna matriz A = (a, y ), 
i,j = 1,...,3, isto é, n = 3 e, em seguida, generalizamos os resultados para qual- 
quer dimensào corno segue: 

Considere LU = A: 


1 

0 

0~ 



U 12 

U 13 


1 

a i2 

a i3 

f, 21 

1 

0 


0 

u 22 

U 23 

= 

a 21 

a 22 

a 23 

co 

_1 

^32 

1 


0 

0 

U 33_ 


_ a 31 

a 32 

a 33_ 


L U A 


Dessa forma, temos urna igualdade envolvendo 0 produto de duas matrizes e 
determinamos as linhas da matriz U e as colunas da matriz L da seguinte forma: 

• 1 3 linha de U 
1 u u = a„ = a„ 

1 u 12 = a 12 —> u 12 = a 12 

1 u 13 = a 13 —» u 13 = a 13 










1 a coluna de L 


fu - 1 


^ 21 U 11 _ a 21 ^ 21 


C 31 U 11 ~ a 31 ^ 31 



/Un 

/Un 


• 2 a linha de U 

-é* 21 ^12 ^22 — ^22 ^ ^22 — ^22 ^ 21^12 
^21 a i3 + ^23 = a 23 ^23 = a 23 ~ ^ 21^13 


• 2 a coluna de L 

( 22 = 1 

/31 ^12 + ^ 32U22 = a 32 ^32 = (&32 ~ ^31^12) ^22 

• 3 a linha de U 

^ 31^13 + 1 32 a 23 + a 33 = a 33 a 33 = a 33 ~ ^ 31^13 ~ ^32^23 

• 3 a coluna de L 

^33 = 1 

Podemos generalizar esse procedimento para urna matriz A = {a ij ), 
i,j = 1,...,n e obtemos: 

Matriz triangular U 

i-i 

U a= a u-^L 'ik u kj U = U.,n i<j. 

k=1 




Matriz triangular L 


ì-i 

a ij ~ S ^ik U kj 

('a = -—- i,j = / > j. 

u n 

Assim, podemos escrever o Algoritmo de Decomposigào, A = LU, corno segue: 

Algoritmo 2.3 

Para m = 1,...,n-1, faga 

Para j = m, m + 1,..., n faga 

m-1 

u m j=a mJ -^( mk u kJ 

k=1 


Para / = m + 1,...,n, faga 


m-1 


a im ^ ik a, 


km 


f. 

•• 


i =1 

1 mm ' 


k=1 


n-1 

ann ~ a nn — ^ (- nk £V, 

k=1 


J kn 


= 1 


Exemplo 2.6 

Considere a matriz A dada por: 


A = 


2 

-4 

0 


1 0 
2 1 
4 2 


Como A ; = Ia„ I = 2*0, A 2 = 



a i2 


2 

1 

a 21 

a 22 


-4 

2 


= 8*0, podemos escrever 


A = LU de forma ùnica e seguindo os passos do Algoritmo 2.3: 












1 a linha de U 


• 1 a colunadeL 

Cu = 1 

1 21 = a 21 / u u = -4 / 2 = -2 
£ 31 = ^31 / Un = 0 / 2 = 0 

• 2 a linha de U 

u 22 = a 22 -1 21 u 12 = (2) - (2)(2) = 4 

U23 = &23 ~ ^21 U-/3 = (1) ~ ( — 2 )( 0 ) = 1 

• 2 a coluna de L 

P 22 = 1 

i 32 = (a 32 - i 31 u 12 )/u 22 = ((4) - mi)) / 4 = 1 

• 3 a linha de U 

U 33 = a33 ~ ^31 U 13 — 132 U 23 = ( 2 ) — mo) — (1)(1) = 1 

• 3 a coluna de L 

^33 = 1 

Assim, obtemos as matrizes triangular interior L e superior U 



1 

0 

0 ' 


~2 

1 

0 ~ 

L = 

-2 

1 

0 

u = 

0 

4 

1 


0 

1 

1 


0 

0 

1 








Aplicacào na resolugào de sistemas de equagòes lineares 


Seja Ax = b um sistema de equagòes lineares, em que A = (a /y ), i,j = 1,...,n 
com det(A) * 0. 

Considerando que a matriz A do sistema satisfaz às condigòes do Teorema 
2.1, temos que A = LU e, portanto, o sistema Ax = b pode ser escrito: 

( LU)x = b 

Nesse caso, denominando Ux = y, a solugào do sistema dado consiste na 
resolugào de dois sistemas triangulares (interior) e (superior), respectivamente, 
corno segue: 

r Ly =b 
I ux = y 


Para resolver os sistemas de equagòes triangulares obtidos, basta usar os 
Algoritmos 2.1 e 2.2, respectivamente. 

Exemplo 2.7 

Considere o seguinte sistema de equagòes lineares: 


2 

1 0 ' 


V 


3 

-4 

0 

2 1 

4 2 


x 2 

_ X 3 _ 

~ 

1 

1 

4^ f\D 

i_ 


Para aplicagào do Mètodo de Decomposigào LU, verificamos as condigòes 
de existència e unicidade da matriz 

A = LU , utilizando o Teorema 2.1 : 

1. Temos que: 


A ,= 2=2*0 


A 2 - 


2 

-4 


1 

2 


= 8*0 












det(A) = 


2 1 0 
4 2 1 


= 8*0 


0 4 2 

Portanto, a matriz A satisfaz condigòes do Teorema 2.1 e det(A) * 0 e po- 
demos resolver o sistema dado usando o Mètodo de Decomposigào LU. 

2. Construpào das matrizes L e U: 

Usando o Algoritmo 2.3 temos a decomposipào da matriz A = LU da se¬ 
guirne forma: 



1 

0 

0 ~ 


~2 

1 

0 ~ 

L = 

-2 

1 

0 

e U = 

0 

4 

1 


0 

1 

1 


_0 

0 

1 _ 


3. Solugào dos sistemas de equagòes triangulares 
a) Ly = b sistema triangular interior 


1 

0 

0~ 

V 


3~ 

-2 

1 

0 


= 

-2 

0 

1 

1 

y 3 _ 


4 


Resolvendo o sistema de equagòes triangular superior obtido, usando o 
Algoritmo 2.1 temos: 

y 1 =3 

< -2y 1 + 1y 2 = 2 y 2 = 4 

y 2 +y 3 = 4 ^y 3 =° 


Assim temos y = (3, 4, 0) { 

b) Ux = y sistema triangular superior 


" 2 

1 

0 


x i 


3 


x i 


1 

0 

4 

1 


x 2 

= 

4 

— > 

x 2 

= 

1 

0 

0 

1 


x 3 


0 


x s 


0 
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Resolvendo o sistema de equagòes triangular superior, usando o Algoritmo 
2.2 temos: 

x 3 =0 

< 4x 2 + x 3 = 4 — > x 2 = 1 
2x 1 + x 2 = 3 —> x 1 =1 

Portanto, a solugào do sistema é dada por x = (1,1,0)’. 

Métodos de Eliminagào 

Os métodos de eliminagào para resolugào de um sistema de equagòes 
lineares Ax = b , consistem em transformà-lo num sistema equivalente triangular 
superior aplicando operagòes elementares sobre linhas corno multiplicar urna 
linha por constante diferente de zero e subtrair de outra linha, as quais serào 
exibidas no mètodo a seguir. 


2.2.4 Mètodo de Eliminagào de Gauss 


Considere o sistema de equagòes lineares Ax = b, em que A = (a (/ ), 
i,j = 1,...,n é matriz dos coeficientes, com det(A) * 0, x = (x p x 2 ,..., xj é o ve- 

tor incògnita transposto e b = (£>,, b 2 ,...,b n ) t é o vetor dos termos independentes 
transposto. 

Representamos o sistema de equagòes lineares da seguinte forma: 


a„x 7 + a 12 x 2 + ... + a 1n x n = b , 
a 21 x 1 + a 22 x 2 + ... + a 2n x n = b 2 


a n1 X 1 + a n2 X 2 + 


+ a nn X n=b n 


O Mètodo de Eliminagào de Gauss, com pivotamento sobre os elementos 
da diagonal da matriz A, consiste em transformar o sistema dado por meio de 
operagòes elementares sobre as linhas, em um sistema equivalente triangular 
superior, tornando, em cada passo, corno pivò, os elementos diferentes de zero 
da diagonal da matriz A. 


(Ab) 


■> (A (n - ,) ,b (n - ?) ) 







em que A (n 1) x = b {n 1] é um sistema triangular superior depois de aplicados 
n-1 passos. 

Consideremos o sistema dado, escrito na seguinte forma: 




,(f). 






a i1 X 1 + a i2 X 2 + a i3 X 3 + a i4 X 4 + ■■■ + a in X n ~ a in+1 










a 21 X 1 + a 22 X 2 + a 23 X 3 + 3 24 X 4 + ••• + a 2r) X n ~ S 2n+1 


~U) V , ~d) v , _(f) v 
d n1 x 1 + d n2 x 2 + d n3 X 3 


, a (1) Y _ -(*) 

• • ■ ‘ d nrv / ' n d nn+1 


Considere a matriz aumentada, isto é, acoplamos o vetor b dos termos 
independentes juntamente com a matriz A do sistema, para que este também 
sofra simultaneamente todas as operagòes de pivotamento, de forma a obter no 
final dessas operagòes um sistema triangular superior equivalente ao sistema 
dado inicialmente. 

Assim, podemos escrever a matriz aumentada da seguinte forma: 



\a (1) 

a ii 

a {1) 

d 12 

a {1) 

d 1n 

a {1) 

d 1n+1 

(A 6) = 

a [1) 

d 21 

a (1) 

d 22 ’ • ' 

a {1) 

d 2n 

a U) 

d 2n+1 


a {1) 

L d n1 

a {1) 

d n2 ■ ■■ 

a (1) 

d nn 

a {1) 

d nn+1 

em que 

a {1) - 

a H ~ 

ti = 1,- 

n+1 

p a {1) 
u d in+1 


Para facilidade de entendimento no procedimento de aplicagòes de opera¬ 
gòes elementares no sistema de equagòes, também chamado de processo de 
pivotamento sobre linhas, consideramos um sistema de equagòes com dimen- 
sào n = 3 e, finalmente, generalizamos operagòes para um sistema de equa¬ 
gòes de qualquer dimensào. 

Considere a matriz aumentada de um sistema de equagòes lineares Ax = b 
com dimensào n = 3 : 







Assim, podemos escrever: 



[g {1) 

a ii 

a {1) 

d 12 

a * (1) 

d 13 

■ a (1) ~ 

d 14 

(Ai>) = 

a (1) 

d 21 

a (1) 

d 22 

a (1) 

d 23 

■ a (1) 

a 24 


a (1) 

a 31 

a {1) 

d 32 

a (1) 

d 33 

■ a 34_ 

em que 

a O) _ 
a ij - 

-- a u , i,J 

= 1,.. 

■ Aeà-4 


Passos do processo de eliminagào 

• Passo 1 

Nesse passo, considere o elemento pivò * 0. Caso este seja igual a 
zero, troca-se as linhas na matriz aumentada, de forma a obter um elemento 
diferente de zero corno pivò. 

Considere os multiplicadores do 1 a passo: 

a (1) 

mii =J k’i = 2’3 

a \ì 


( 0 ) 0 )\ 

e as operagòes de pivotamento sobre as linhas de [A, b) do tipo: 


gO) 

a)Multiplica-se a 1 a linha fixa da matriz aumentada pela constante —e 

a ii 

subtrai da 2 a linha; faz-se a substituigào na 2 a linha por està modificada 
dessas operagòes. 


b) Multiplica-se a 1 a linha pela constante 


. o ) 

*31 

OD 

*11 


e subtrai da 3 a linha; faz-se a 


substituigào na 3 a linha por està modificada dessas operagòes. 


Dessa forma, conservamos a 1 a linha e tornamos nulos todos os elementos 

da 1 a coluna abaixo da diagonal na matriz aumentada (A m , b {1) ), pois ao efetuar 
essas operagòes de pivotamento temos: 




~ 0 ) _ a 21 
d 21 


aO) _ n 

( 1 ) a H ~ U 


a O) _ a 31 
d 31 


a {1) -n 

(1) a H ~ U 


l 11 


1 11 










Efetuando o mesmo càlculo para toda a linha, obtemos os elementos modi- 
ficados, os quais sào marcados corno a‘ y y) , e realizamos o mesmo procedimento 
para os elementos das linhas abaixo da diagonal: 

Assim, escrevemos a matriz aumentada modificada corno segue: 


{A a \b (1) ) 


72 

a {1) ■ 

“73 

a [1) 

d 14 

» ( 2) 

a (2) - 

a (2) 

‘22 

a 23 

a 24 

,(3) 



l 32 

a 33 

a 34 


De forma geral, para um sistema Ax = b , em que A = (a /y ), i,j = te- 

a (o 

mos a (2) = a] 1 / - m j1 a ( ,J\ m n = / = 2,...,n, j = 1,2,..., n + 1. 

a u 

• Passo 2 

Considere agora o coeficiente a (2 J * 0 corno elemento pivò. Caso = 0, 
troca-se com as linhas seguintes até que o coeficiente afj que ocupa a segun- 
da linha e segunda coluna seja diferente de zero. 

Considere os multiplicadores do 2 2 passo: 


p (2) 

a 12 


e faga as operagòes sobre as linhas na nova matriz aumentada ( A (1) , b (1) ) modifi¬ 
cada com as operagòes de pivotamento do 1° passo: 

a (2) 

a) Multiplica-se a 2 a linha fixa por e subtrai da 3 a linha; faz-se a subs- 

d 22 

tituigào na 3 a linha por està modificada dessas operagòes. 

Dessa forma, tornamos nulos todos os elementos da 2 a coluna abaixo da dia¬ 
gonal na matriz (A {2 \ b {2) ), pois ao efetuar essas operagòes nessa matriz temos: 


a 


( 2 ) 

32 


,( 2 ) 

*32 (2) 


,( 2 ) d 22 

122 


= 0 


Efetuando o mesmo càlculo para toda a linha, temos os elementos modifi- 
cados, os quais serào marcados corno af\ e procedemos igualmente para as 
linhas abaixo da diagonal: 
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Assim, podemos escrever a matriz modificada da seguinte forma: 



1 

a {1) 

d 12 

a {1) 

d 13 

a (1) 

d 14 

II 

0 

P (2) 

d 22 

a (2) 

d 23 

a (2) 

d 24 


0 

0 

a {3) 

d 33 

a {3) 

d 34 


De forma geral, para um sistema Ax = b, em que A=(a iJ ), i,j = 1,...,n, 
temos: 


o(3> 

d i j ~ d i j ' n i2 d 2j 


( 2 ) 


m, 


i (2) 

3 l2) 

a 22 


3,...,n j = 2 ,...,n + 1 


Assim, depois de executados 2 passos, obtemos o sistema inicial dado 
Ax = b na forma equivalente triangular superior, da seguinte maneira: 




^b 


,(f). 


a i1 X 1 + a i2 X 2 + a i3 X 3 


,( 2 ), 


,(2k 


- a {1) 

- a 14 

,( 2 ) 


a 22 X 2 + a 23 X 3 — a 24 


,( 3 ), 


,( 3 ) 


a 33 X 3 ~ a 34 


A solupào é dada, conforme anteriormente, usando a expressào geral para 
a solupào de um sistema triangular superior: 


X: = 


i = i+i 


a ij X j 


a ii 


i = 3, 


1 


Generalizando essas operagòes para um sistema com n linhas e n colunas, 
necessitamos de n-1 passos para transformà-lo em sistema equivalente trian¬ 
gular superior. Essas operagòes sào detalhadas no Algoritmo de Eliminagào de 
Gauss, descrito corno segue: 









Algoritmo 2.4 


a)Construgào do sistema triangular superior equivalente 
Para k = faga 

Para / = k + 1,...,n, faga 


m 


a (k) 
m = _ 

ik oW 
d kk 


Para j = k,..., n + 1 


+ 1) _ q(^) 


U 


'1 


_ (k) (k) 

'"ik d kj 


b)Calcular a solugào do sistema triangular superior 

Usar o Algoritmo 2.2 para resolugào de sistemas de equagòes triangular 
superior. 

Exemplo 2.8 

Considere o seguirne sistema de equagòes lineares: 


"3 

2 

1 

~ x i~ 


"4" 

1 

1 

1 

x 2 

= 

1 

2 

1 

1 

. X 3_ 


2 


Como det (A) = 


3 

1 

2 


2 

1 

1 


1 

1 

1 


= 1*0, o sistema possui urna ùnica solugào, 


usando o Mètodo de Eliminagào de Gauss, temos: 

Considere a matriz aumentada, conforme a seguir: 


(A b) = 


3 2 1 
1 1 1 
2 1 1 


4 

1 

2 


Executando o Passo 1 e o Passo 2, respectivamente, temos: 

• Passo 1 

Considerando que o pivò afl = 3*0, podemos calcular os multiplicadores: 
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m„ = 


,( 7 ) 

*21 

,(?) 

*11 


1 _ 

3 


e m 


21 


a (1) 

_ d 31 _ 

.(») 

*11 


2 _ 

3 


Multiplica-se a 1 a linha por 1/3 e subtrai da 2 a linha; faz-se a substituigào 
na 2 a linha, em seguida multiplica-se a 1 a linha por 2/3 e subtrai da 3 a linha e, 
novamente, faz-se a substituigào na 3 a linha. 

Assim temos a matriz aumentada modificada da seguinte forma: 


(A { 1 \b {1) ) = 


3 

0 

0 


2 1 
1/3 2/3 

-1/3 1/3 


4 

-1/3 

-2/3 


• Passo 2 

Como o pivò a 22 = 1/3 * 0, calculamos os multiplicadores corno segue: 


m 32 ~ 


,( 2 ) 
d 32 




-1/3 

1/3 


= -1 


Multiplica-se a 2 a linha por -/, subtrai da 3 a linha e faz-se a substituigào na 
3 a linha. Temos a matriz aumentada modificada corno segue: 


(A { 2 \b (2) ) = 


3 

0 

0 


2 1 

1/3 2/3 
0 1 


4 

-1/3 

-1 


Assim, obtemos um sistema de equagòes triangular superior, o qual pode 
ser reescrito da seguinte forma: 

3x 1 + 2x 2 +1x 3 =4 

1/3x 2 + 2/3x 3 = - 1/3 
1x 3 = -1 


Resolvendo o sistema triangular superior equivalente, usando o Algoritmo 
2.2, temos: 

1x 3 — — 1 —^ x 3 = — 1 

< 1/3x 2 + 2/3x 3 = -1/3 -» x 2 = 1 
3x 1 + 2x 2 +1x 3 =4 x 1 = 1 














Portanto, a solugào do sistema é dada por x = ( 1,1, -1) { . 

Observagào 

Para verificar-se que a solugào obtida é solugào do sistema originai dado, 
basta substitui-la nas equagòes do sistema e verificar a igualdade. 

Aplicagào 

Urna aplicagào de resolugào de sistemas de equagòes lineares é quando 
calculamos a matriz inversa de urna matriz A, em que A = (a jj ), i,j = 1,... ,n. Nesse 
caso necessitamos da resolugào de n sistemas de equagòes lineares, em que 
usamos os métodos diretos vistos anteriormente. 

2.3 Inversas de matrizes 


Considere A = (a ,,), i,j = 1,... ,n urna matriz nào singular, isto é, det(A) * 0. 
Entào existe urna ùnica matriz A 1 chamada de inversa de A, de tal modo que 
AA _1 =1, em que a matriz I é urna matriz diagonal, conhecida corno a matriz 
identidade, cujos elementos da diagonal sào todos iguais ao nùmero 1. 

Dessa forma temos: 


f *\1 

a 12 . . 

• 

( X 

A 11 

X 12 ■ 

■ X 1n 


fi 

0 

0 . 

. 0 ' 

a 21 

a 22 ■ ■ 

• a 2n 

X 21 

X 22 ■ 

■ X 2 n 

= 

0 

1 

0 

. 0 

y^nl 

a n2 ■■ 

^nn j 

K X n1 

X n2 •• 

• X nn y 


° 

0 

0 . 

. 1 

J 


Portanto, para determinar as n colunas da matriz inversa A 1 , temos que 
resolver n sistemas de equagòes lineares, usando qualquer um dos métodos 
diretos vistos anteriormente. 

Para calcular a 1 a coluna da matriz inversa A -1 , resolvemos o seguirne sis¬ 
tema de equagòes lineares: 



a 


12 


a 


21 


a 


22 



a 


n2 


In 


2n 


'X ' 
11 


21 


^ 

0 


o 

v y 
















Assim sucessivamente, calculamos a n-ésima coluna da matriz A \ resol- 
vendo o seguinte sistema de equagòes lineares: 


a a 

11 12 


a a 

21 22 



(x 

In 


\ 



'o N 

0 


v 


a 


ni 


a 


n2 


a 

nn 


X 

\ nn J 


y 


Exemplo 2.9 

Considere a matriz A dada por: 


A = 


1 

1 

0 


0 2 
1 0 
0 2 


Temos que a matriz inversa de A pode ser escrita corno: 


X 

X 

X 

11 

12 

13 

X 

X 

X 

21 

22 

23 

X 

X 

X 

31 

32 

33. 


Usando o Mètodo de Eliminagào de Gauss, podemos calcular a inversa da 
matriz A dada, corno segue: 


Como det(À) 


1 0 2 
1 1 0 
0 0 2 


= 2*0 temos que a matriz inversa existe e é ùnica. 


Considerando que AA 1 =1, podemos escrever: 


1 

0 

2 " 


1 

1 

0 


0 

0 

2 



xxx 

11 12 13 

XXX 

21 22 23 

XXX 

31 32 33. 


1- 

1_ 

0 

1 - 

O 

= 

0 

1 

0 

-1 

1 

O 

0 

1 


Temos très sistemas de equagòes lineares para serem resolvidos, corno 
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seguem: 






















1 a coluna da matriz inversa de A, 


1 0 2 


X 


V 

1 1 0 


ii 

X 


0 


21 


0 0 2 


1 

>< 

Co 

1_ 


0 


2 a coluna da matriz inversa de A, 


~1 0 2 


X 


~0~ 



12 



1 1 0 


X 

— 

1 



22 



0 0 2 


X 

L 32 J 


0 


3 a coluna da matriz inversa de A, 


1 0 2 


X 


" 0 " 



13 



1 1 0 


X 

— 

0 



23 



0 0 2 


X 

L 33 J 


1 


Construimos inicialmente a matriz aumentada (A, I) formada pela matriz A 
e os termos independentes dos sistemas obtidos anteriormente. Aplicamos as 
operagòes do Mètodo de Eliminagào de Gauss, transformando a matriz dada na 
forma triangular superior, e os vetores dos termos independentes sào modifica- 
dos, corno segue: 


' 1 

0 

2 

1 

0 

0 ~ 


' 1 

0 

2 

1 

0 

0 ~ 

1 

1 

0 

0 

1 

0 

~ 

0 

1 

-2 

-1 

1 

0 

0 

0 

2 

0 

0 

1 


0 

0 

2 

0 

0 

1 


Assim, podemos resolver os sistemas triangulares obtidos: 


' 1 0 2 


X 


f 



ii 



0 1 -2 


X 

— 

-1 



21 



0 0 2 


1 

CO 

>< 

_ 1 


0 


Resolvendo o sistema de equagòes triangular superior obtido, temos a so- 
lugào da 1 a coluna da matriz inversa de A, dada por: 


(x„ x 21 x 31 i = {1-1,0)' 
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Para determinar a 2 a coluna da matriz inversa de A, resolvemos o sistema 
de equagòes triangular superior obtido: 


~1 0 2 


X 


~ 0 ~ 



12 



0 1 -2 


X 

— 

1 



22 



0 0 2 


1 

l\3 

1_ 


0 


cuja solugào é dada por: 

(Xy 2 X22 X 32 ) = iPì 0) 

Finalmente para determinar a 3 a coluna da matriz inversa de A, resolvemos 
0 sistema de equagòes triangular obtido corno segue: 


1 0 2 


X 


~0~ 



13 



0 1 -2 


X 

— 

0 



23 



0 0 2 


X 

L 33 J 


1 


cuja solugào é dada por: 

{x 13 x 23 x 33 y = (- 1,1,1/2)' 


Portanto, temos a matriz inversa A" 1 da matriz dada por: 


A 1 


1 0 -1 

-1 1 1 

0 0 1/2 


Observagào 

Quando usamos 0 Mètodo de Decomposigào LU no càlculo da matriz inver¬ 
sa, transformamos a matriz A no produto de duas matrizes triangulares, interior 
e superior respectivamente. Nesse caso, teremos très sistemas de equagòes li- 
neares para serem resolvidos e, consequentemente, dois sistemas de equagòes 
triangulares inferiores e superiores em cada um deles. 


















2.4 Métodos iterativos para resolugào de sistemas de equagoes 
lineares 

Podemos resolver um sistema de equagoes lineares Ax = b utilizando para 
isto métodos iterativos, os quais consistem em determinar urna sequència de 
solugòes aproximadas, em que cada solugào obtida depende da solugào ante- 
rior pela aplicagào do mesmo procedimento. 

Para construir um mètodo iterativo, escrevemos o sistema de equagoes 
lineares Ax = b, na forma equivalente x = Hx + g. Essa forma equivalente sera 
detalhada no Mètodo iterativo de Jacobi-Richardson, conforme segue. 

De maneira geral, obtida a forma equivalente x= Hx + g, a partir de urna 
solugào aproximada inicial x (0) , determinamos a sequència de solugòes aproxi¬ 
madas considerando o processo iterativo: 

x U) = Hx (0] +g 
x (2) = Hx (1) +g 


x (k+1) = Hx {k) +g, k = 0,1,2,... 

Nesse caso, temos H = (/? /y ), i,j = 1,..., n chamada de matriz iterativa e um 
vetor g = (g,), i = 1,...,n. 

Assim, a partir de urna solugào aproximada inicial x (0) = {x\°\ xx [ n 0) ) 
para a solugào exata x = (x 1 ,x 2 ,...,x n ) do sistema Ax = b, determina-se urna 

sequència de vetores x {1 \ x {2) , x°\... e espera-se que seja convergente para a 
solugào x, isto è, 

jim x (k) = x 


Apresentamos a seguir alguns resultados e definigòes necessàrios para 
entendimento do Mètodo iterativo de Jacobi, corno seguem: 

Definigào 2.2 Norma de vetores 

A norma de um vetor x e V (espago vetorial) è urna fungào ||. || ; V -» 9t + , 
satisfazendo às seguintes condigòes: 

a) | x II > 0; | x II = 0 <^> x = 0, V x e V 




b) | ax | = | a || x ; Vae 9t, Vxe V 

c) ||x + y |<|x|+|y ||;Vx, ye V 


Exemplos de normas: 

Considere o espago 9?" formado pelos vetores x = (Xp x 2 ,..., x n ): 
Temos as seguintes normas de vetores: 

III I | | n | | 

X II = X + X + ... + X =Y X 
11 / | 1 | | 2 | | n | 4 - J \ i | 




( n „ X 

X 1 =lx 2 + X 2 + . 

Il 2 \ 1 2 

, + x 2 = 

n 

xi -, r 



V’ 1 J 


X 


max ix 

, x L 

.. Ix f 

[ = max x. 

Il 1 1 

1 2 \ 

1 n [ 

J 1<i<n i \ 


Essas normas sào chamadas de norma 1, norma 2 ou Euclidiana e norma 
infinita, respectivamente. 

Observagào 

Podemos interpretar as normas definidas anteriormente no 3ì 2 , bastando, 
para isso, representar o ponto x = (x 1 ,x 2 ) t no plano e observar, por exemplo, que 
a norma 2 é a distància entre o ponto x e a origem do sistema de coordenadas. 

Exemplo 2.10 

Considere x = (x p x 2 , x 3 ) = (-1,2,-4)& 9t 3 . 

Calcolando as normas do vetor x, obtemos: 


x 



-1 + 2 + 


4 


= 7 



( 3 


2 N 

l x L = 


X 

/ 





) 


y/l 2 +2 2 +4 2 




Il II 

, r 

X [ 

o = maxi 



} = max {7, 2, 4} = 4 


Definigào 2.3 Sequència convergente 

Dizemos que urna sequència de vetores x (,) = (x p x 2 ,..., x n ) 0) e 9i n con- 




















verge para o vetor x = (x p x 2 ,..., x n ) e 9T se 


x (/) -x 


->0, quando /—>°°, para qualquer norma em SR n . 


Exemplo 2.11 

Considere a sequència x (,) = 


i 


e9t n , i = 1,2,... 


V V 

Temos que, à medida que /—>°°, a sequència converge para o vetor 
x = (0,0,..., 0). 


Ou ainda, 
Como 


x (/) -x 


= -, temos 


x (/) -x 


-> 0, quando / —^ 00 . 


Definigào 2.4 Norma de matriz 

Considere o espago SR(n, n) das matrizes quadradas A = (a,.) /, y = 7,..., n. 
Definimos as seguintes normas de matrizes: 


I A II =\\A IL = maxY 

i li; il no l<j< n “ 


I /4 II = Il II, = max Y 


1<i<n 


H 


Essas normas sào chamadas norma coluna e norma linha, respectivamente. 

Exemplo 2.12 

Considere a matriz A: 


A = 


7 2 3 
1 1 1 
1 2 4 


Entào, 


A 



max Y 

y </< 3 % 


a 

ij 


màx\3, 5, 8} = 8 

1<j<3 L J 


























>A|| =\\A\\=màx'Y a = màx\6,3,7} = 7 

Ile*, Il IIL 1<i<3 ' Il islsn L J 


H 


1<i<3 


Definicào 2.5 Normas consistentes 

Considere urna norma de um vetor xe9t n e urna norma de urna matriz 
Ae 9?(n, n). Dizemos que essas normas sào consistentes se: 

| Ax | < | A II x ||,Wte 9?(n,n) e Vxe'K". 

Propriedade 

Il /Ae II < Il /A mi e ||, VASe 9t(n,n). 

Definicào 2.6 Matriz estritamente diagonalmente dominante 

Dizemos que urna matriz A = (a /; ), i,j = 1,...,n é estritamente diagonal¬ 
mente dominante se: 



n 


a 

// 

V 

M 

a 

ij 

H. 


'W 


Exemplo 2.13 


A = 


-4 

1 

1 


1 

-2 

1 


0 

0 

3 


Usando a definigào 2.6, temos que: 

/ = /-»|-4| > |y| + |0| 

/ = 2 -» | - 2 | > | /1 +1 0 | 

/ = 3—»|3|>]f| + |f[ 

Portanto, concluimos que a matriz dada é estritamente diagonalmente 
dominante. 












2.4.1 Mètodo iterativo de Jacobi-Richardson 


Considere o sistema de equagòes lineares Ax = b, em que A = {a jj ), 
i,j = é a matriz dos coeficientes com det(A) * 0 e com os elementos da 
diagonal principal a„ *0, i = 1,...,n, x = {x 1 ,x 2 ,...,xj é o vetor incògnita trans¬ 
posto e b = (b p b 2 ,... ,b n )’ é o vetor dos termos independentes transposto. 

Representamos o sistema dado na seguirne forma: 

a n x i + a i 2 x 2 + a 13 x 3 + ... + a 1n x n = b , 
a 2 i x i + & 22 x 2 + a 23 x 3 + ■ • ■ + 3 2n x n = b 2 

< 

a n1 X 1 + a n2 x 2 + a n3 X 3 + - + a nn X n = b n 

O Mètodo iterativo de Jacobi consiste em escrevermos o sistema Ax = b 
na forma equivalente x = Hx + g, corno segue: 

Dividindo cada linha pelo elemento da diagonal e explicitando x 1 na 1 a equa- 
gào, x 2 na 2 a equagào, x 3 na 3 a equagào e até x n na n-ésima equagào, temos: 

X 1 = — ( b 1 - a i2 X 2 ~ a i3 X 3 ~ ■ • ■ — a in X n ) 

a i1 

X 2 = ( b 2 ~ a 21 X 1 ~ a 23 X 3 ~ ■ ■ ■ ~ a 2n X n ) 

a 22 

X n = ^~ (b n -a n1 x 1 - a n2 x 2 - ...-a nn _ f x„_,) 

l a nn 

Na forma matricial temos: 


V 


f 

0 

a 21 

a 31 

_ a n1 

\ 






a i1 

a i1 

a i1 




a u 

X 2 

= 

a 21 

a 22 

0 

a 23 

a 22 

a 2n 

a 22 


X 2 

+ 

È2_ 

a 22 



_ a n1 

a n2 


0 




bn_ 



v a nn 

a nn 


/ 



y a nn À 














Assim, podemos escrever x = Hx + g, e m que 


h a=< 


0, i = j 


a ij ■ ■ 

i*J 

a, 7 


é chamada de matriz iterativa e g j = —, / = 

dii 


Podemos escrever o processo iterativo de Jacobi corno segue 
x (k+1) = Hx ik) +g k = 1,2,... 

Assim: 


f xS k+1 i \ 


x 2 ( k+1 ) 


x (k+r) 
v n y 


0 

a 27 

a 31 

_ a n1 

a 7 7 

a 77 

a i 1 

a 27 

0 

a 23 

a 2r> 

a 22 

a 22 

a 22 

fnl 

a n2 


.. 0 

a nn 

a nn 



! K A 


'x/*) A 


x 2 (' < ) 


X„(*> 

\ n J 


bi 


*22 


v a "V 


Podemos, ainda, escrevè-lo na seguirne forma: 

X/k+D = 


(*)_ ^X/) 
a 7? a 7/ 


. - j*) + A. 


-*7 7 


-^7 7 


x 2 ( ~ k+1 '> = 


- ^Lx/'O - ^X 3 (^) 
a 22 a 22 


X (k+o = - ^2Lx/*) - —XJ k 1 - 


hn_x W + A_ 

n 

a 22 a 22 


- ^2=Lx n ,(*)+ Al 

a n/7 " a nn 















Estudo da convergència 

Teorema 2.2 

Seja urna norma matricial consistente com alguma norma vetorial e 
x (0) e9T urna solugào aproximada inicial qualquer para o sistema de equa- 
gòes. Se | H \\<1, entào a sequència de solugòes aproximadas definida pelo 
processo iterativo x [k+1) = Hx {k) +g, k = 0, 1 , 2 ,... converge para a solugào x 
do sistema de equagòes Ax = b. 

Prova: Referència [1]. 

Resultado: 

É possivel provar que se a matriz A do sistema Ax = b for estritamente 
diagonalmente dominante, entào o Mètodo de Jacobi-Richardson gera urna se¬ 
quència de solugòes aproximadas convergente para a solugào do sistema, pois 
sempre temos || H ^ < 1. 

Apresentamos de forma detalhada o Algoritmo do Mètodo de Jacobi-Ri¬ 
chardson, observando que o critèrio de parada na sequència gerada de vetores 
è o erro relativo visto anteriormente, apenas usamos a notagào de norma, pois 
estamos trabalhando com sequència de vetores. 

Algoritmo 2.5 

a) Fornega urna solugào inicial aproximada x (0) = ( x\ 0) , x { 2 0) ,..., x ( n 0) ) e e > 0 
urna tolerància fixa. 

Faga k = 0 e PARE = FALSO 

b) Construgào da sequència de solugòes aproximadas: 

Enquanto PARE = FALSO, faga: 

Para i = faga: 


x 


(k+1) 


~ a ij X T + I - a a x 

i=i j=i+i 


\ 


( k ) 
i 


J 


/a ii + y 
a ii 


X (k+1 ) 

Se -il —tt-tt—^< e, entào PARE = VERDADE 

X (k+1) 


Senào k = k + 1 
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Exemplo 2.14 

Considere o seguinte sistema de equagòes lineares: 


~4 

1 

1~ 

x i 


" 3" 

1 

4 

1 

x 2 

= 

0 

_1 

1 

5_ 

_ X 3 _ 


-4 


Usando o Mètodo iterativo de Jacobi, podemos a partir de urna solupào 
aproximada inicial x (0) = (x y (0) , x { °\ x ( 3 0) )' =(0,0,0)' calcular urna solupào apro- 
ximada para o sistema dado com urna precisào e = 10'. 

Assim, a construgào da matriz iterativa H é dada por: 


H = 


0 

1/4 

1/5 


-1/4 

0 

-1/5 


-1/4 

-1/4 

0 


e para verificar a condipào de convergència temos: 


H II = Il H II = max VI /?, I = màx{2/4, 2/4, 2/5} = 0.5<1 

\\L II lloo i \ J | j L J 


Portanto, a sequència de solupòes aproximadas é convergente para a so- 
lupào do sistema dado. 

Càlculo das iteragòes: 


^(k+1ì _ 




y(k + 1) _ 

- 


X 


(k+1) _ 
3 


X 


( k) 


— X 


( k) 


— X 


( k ) 


— X 

4 


( k) 


— X 

4 


(.k) 


— X 


( k ) 


3 

+ — 

4 

+ 0 

4 

5 


k = 0,1,2,... 


Para k = 0 e tornando x (0) = (0,0,0)', temos x 


|x (? >-x<°> 
Il y-m II 


(?) _ 


V 


-, 0 ,-~ 

4 5 
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Teste de parada: 


= 1>10 
















Para k = 1, temos: 


x {2) = 0.9500 
xf = 0.0125 
xf = -0.9500 

Portanto, 

x< 2 > = ( 0.9500, 0.0125,-0.9500Y 

Il x < 2 ) -x {1) Il 

n-jj—^ = 0.2000 > IO- 1 

x< 2 > 


x\ 2) = 


A2) _ 


A 2 ) - 


1 ( 7 ) 1 ( 1 ) ^ 

—x‘ 7) - x (1) + — 

4 2 4 3 4 

--x f (,) - -x<’> + o 

4 1 4 3 

- —x {1) - -x (1) - — 

5 1 5 2 5 


—> 


—> 


—» 


Assim, sucessivamente, calculamos: 
Para k = 3, temos: 


x ? (3) = 


y(3) _ 
A 2 


y(3)_ 

A 3 


~-x (2) ~-X (2) +- 

4 2 4 3 4 

--x\ 2) --x {2) +0 

4 1 4 3 

-L X (2)_L X (2)_± 

5 1 5 2 5 


—> 


—> 


—> 


x ? (3) = 0.9844 
xf = 0.0000 
xf =-0.9925 


Portanto, 

x (3) =(0.9844, 0.0000,-0.9925)’ 

| X (3) _x (2) | 

-—tt-tt —— = 0.0400 < IO- 1 , entào pare. 

x (3) 


Logo, x = x (3) = (0.9844,0.0000,-0.9925 ) f é a solupào aproximada para 
sistema dado com a precisào e = 10 








Observagào 

Podemos observar que a sequència de solugòes aproximadas 
x (1) , x (2) , x (3) , x (4 \... converge para a solugào do sistema de equagòes 
x = {1,0,-1) t . 

Observagào 

Os métodos iterativos devem ser aplicados na resolugào de um sistema de 
equagòes lineares Ax = b quando este apresentar a matriz A de grande porte e es- 
parsa, isto é, urna matriz com muitos elementos em que a maioria deles sào nulos. 

Nesse caso, quando determinamos a matriz iterativa H, està contém muitos 
elementos nulos, o que diminuì o esforgo computacional. Caso contràrio, deve- 
mos utilizar os métodos diretos vistos anteriormente. 


2.5 Exercfcios 

1. Considere o seguinte sistema de equagòes lineares: 


'-5 

1 

1 " 

V 


2 

0 

1 

0 

x 2 

= 

1 

2 

3 

2 

. X 3. 


5 


a) Verifique as condigòes de aplicabilidade do Mètodo de Decomposigào 
LU para o sistema dado. 

b) Caso possivel, resolvè-lo usando o Mètodo de Decomposigào LU. 

2. Determinar a inversa da matriz A do exercicio 1), usando o Mètodo de 
Decomposigào LU. 

3. Resolver o seguinte sistema de equagòes lineares usando o Mètodo de 
Eliminagào de Gauss. 


~0 

6 

-1 

2 

V 


f 

0 

5 

-1 

1 

x 2 


0 

2 

1 

-1 

-1 

*3 


-2 

0 

-1 

-2 

1 

_*4. 


-1 


4. Desenvolva o Mètodo de Gauss-Jordan para resolugào de sistemas de 
equagòes lineares e faga um exemplo. Referència [1]. 
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5. Considere o sistema de equagòes lineares: 


~0 

1 

1~ 

V 


~0~ 

2 

3 

0 

x 2 

= 

3 

3 

1 

1 

_ X 3_ 


0 


a) Faga trocas de linhas e, caso possivel, resolva o sistema usando o Mè¬ 
todo de Decomposigào LU. 

b) Determinar a inversa de A, usando o Mètodo de Gauss-Jordan. 

c) Faga trocas de linhas e, em caso de convergència garantida, resolva o 
sistema usando o Mètodo Iterativo de Jacobi com urna solugào inicial 
dada e e = 0.1. 


6. Considere o seguinte sistema de equagòes lineares: 


-9 

2 

2 

V 


"-5" 

1 

3 

0 

x 2 

= 

4 

-1 

1 

3 

X 3. 


3 


Caso haja convergència garantida para o Mètodo Iterativo de Jacobi, resol¬ 
va o sistema dado a partir da solugào inicial x(°) = (0 00) e e = 0.1. 

7. Mostre que || x || m < || x || < n|| x ^, x e R n . 

8. Considere duas cidades C 1 e C 2 com produgòes de lixo de 8 e 5 tonela- 
das, respectivamente. Numa àrea apropriada existem dois aterros A 1 e A 2 com 
capacidades de processar 7 e 6 toneladas de lixo, respectivamente. Deseja-se 
saber quais as quantidades de lixo devem ser enviadas das cidades C 1 e C 2 para 
os aterros A 1 e A 2 de forma que todo o lixo seja transportado. 

a) Escreva a formulagào matemàtica correspondente do problema dado. 

b) Fìesolvè-lo usando o Mètodo de Eliminagào de Gauss. 

c) Interprete os resultados obtidos. 

d) Resolva o problema dado, com o Software Matlab, usando o Mètodo de 
Gauss-Jordan. 

e) Resolva o problema dado, com o Software Numèrico, usando o Mètodo 
de Eliminagào de Gauss - Referència [1]. 
















9. Faga um mapa conceitual detalhado sobre solugào numèrica de sistemas 
de equagòes lineares, considerando os métodos vistos, introduzindo labels indi¬ 
cando se aprendeu (A), nào aprendeu (N), se gostou (G), detestou (D), achou 
interessante (I), etc. 




UNIDADE 3 


Solugào Numèrica de Equagòes 



3.1 Introdugào 


Nas diversas àreas cientificas, frequentemente deparamo-nos com problemas 
reais envolvendo a resolugào de equagòes, isto é, dada a equagào f(x) = 0, dese- 
jamos determinar a solugào (ou solugòes) reai x , tal que f(x) = 0. Por exemplo, 
consideremos a equagào f(x) = sen(x) - ln(x) + 2 = 0, desejamos determinar a 
solugào reai x, tal que f(x) = sen(x)-ln(x) + 2 = 0. 

Métodos iterativos sào desenvolvidos para determinar aproximadamente 
essa solugào reai x, embora tenhamos métodos iterativos especificos para de¬ 
terminar a solugào x, quando està é um nùmero complexo. 

Métodos iterativos sào apresentados para determinar a solugào x, quan¬ 
do està é um valor reai, e para isso, necessitamos de urna solugào aproximada 
inicial x 0 , a partir da qual geramos urna sequència {x n }, n = 1,2,... de solugòes 
aproximadas, que sob determinadas condigòes teóricas converge para a solu¬ 
gào desejada x. 

A solugào inicial x 0 pode ser obtida através de recursos gràficos, no qual 
localizamos urna vizinhanga ou um intervalo [a,b] em que encontra-se a solu¬ 
gào x, conforme exibimos na Figura 3.1 : 



Observando a Figura 3.1, vemos que a solugào x tal que f(x) = 0 encon¬ 
tra-se onde a fungào f(x) corta o eixo das abscissas, isto é, um intervalo em 
que os valores da fungào possui sinais opostos. 

Podemos, entào, tornar urna solugào inicial x 0 nas vizinhangas dessa raiz, 
isto é, no intervalo [a, b\ para inicializar a sequència de solugòes aproximadas du¬ 
rante aplicagào dos métodos iterativos que serào apresentados posteriormente. 









Para ilustrar a solugào numèrica de equagòes, considere o problema de 
nivel de oxigènio em um rio, conforme Figura 3.2 



Figura 3.2 Rio contaminado. 

Sabendo-se que para calcular o nivel de oxigènio N (mg/L) em um rio a 
jusante de urna carga de esgoto, podemos usar a seguinte fungào: 

N(x) = 10-20{e-°- 15x -e-° 5x ) 

em que x è a distància a jusante em quilòmetros, conforme Figura 3.3: 



Figura 3.3 Nivel de oxigènio de um rio. 

a) Deseja-se determinar a jusante em que o nivel do oxigènio cai pela primei- 
ra vez até urna leitura de 5 mg/L, observando-se quais niveis de oxigènio 
abaixo de 5 mg/L sào prejudiciais aos peixes e ao ecossistema aquàtico. 

b) Deseja-se determinar a distància a jusante na qual o oxigènio està no 
valor minimo. Qual a concentragào de oxigènio nessa posigào? 









Definicào 3.1 

Dizemos que x é urna raiz, ou um zero, da fungào f(x) se f(x) = 0. 

Exemplo 3.1 

Seja f(x) = x 2 - 5 = 0. Termos que as raizes da equagào x 2 - 5 = 0 sào: 
x = ± yfiT = 2.236067977 , e nesse caso f(x) = 0. 

3.2 Localizagào das raizes: Mètodo Grafico 

Para localizar urna vizinhanga para a raiz de f(x), tragamos o gràfico de f(x); 
assim, os pontos em que este corta o eixo das abscissas nos informam a raiz (as 
raizes) de f(x). 

Considere o exemplo anterior, f(x) = x 2 -5, e temos as raizes x = ±45, 
conforme gràfico ilustrado na Figura 3.4: 



Figura 3.4 Raizes de urna fungào. 

Observando a Figura 3.4, vemos que o gràfico de f(x) permite identificar 
onde estào aproximadamente as raizes de f(x). Nesse caso, temos as raizes 
x = ± 45 = 2.236067977 com f(x) = 0. 

Podemos ainda transformar a equagào f(x) = 0 na forma equivalente, na 
igualdade de outras fungòes /)(x) = f 2 (x). Os pontos de intersegào dos gràficos 
/)(x) e 400 serào as raizes procuradas, conforme exemplo a seguir: 


Exemplo 3.2 

Considere a equagào f(x) = 4x-e x =0. 
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Podemos escrever a equagào dada na forma equivalente 4x = e x , isto é, 
f 1 (x) = f 2 (x), com f 1 (x) = 4x e f 2 (x) = e x , conforme gràfico da Figura 3.5: 



Figura 3.5 Intersecgào de gràficos de fungòes. 

Observando a Figura 3.5, vemos que a raiz x encontra-se na intersecgào 
dos gràficos de /j(x) e f 2 (x). 

3.3 Métodos numéricos para resolugao de equagoes 

Nesta unidade apresentaremos alguns dos principais métodos para resol¬ 
ver numericamente urna equagào. 

3.3.1 Mètodo da Bissegào 


Considere a fungào f : 9? 91 urna fungào continua. Desejamos resolver a 

equagào f(x) = 0, isto é, determinar urna solugào x reai tal que f(x) = 0. 

O Mètodo da Bissegào é baseado no Teorema do Valor Intermediàrio, o qual 
afirma que se urna fungào é continua no intervalo [a,b], e satisfaz a condigào 
f(a)f(b)<0, valores de f(a) e f(b) com sinais opostos, entào existe xe (, a,b ) 
tal que f(x) = 0, isto é, existe ao menos urna raiz no intervalo \a,b], conforme 
ilustrado na Figura 3.6. 







Figura 3.6 Intervalo contendo urna raiz da fungào f(x). 

O Mètodo da Bissegào consiste em localizar inicialmente um intervalo 
[a,b], em que se encontra a raiz x, e determinar urna sequència de intervalos 
[a/ = 0,1,..., em que r 0 =a e s 0 = b, de forma que a amplitude do intervalo 
numa iteragào seja igual à metade da amplitude do intervalo anterior e que sem¬ 
pre contenha a raiz x. 

A sequència de intervalos sera calculada até que a amplitude do intervalo 
seja menor do que urna tolerància s preestabelecida. 

As sequèncias r j; s ; e x, sào construidas a partir de um intervalo inicial 
[r 0 , s 0 ] tal que f(r 0 )f(s 0 ) < 0, conforme Figura 3.7: 



Figura 3.7 Intervalo inicial contendo urna raiz da fungào f(x). 














Determina-se o ponto mèdio do intervalo [r 0 , s 0 ], dado por: 


x 0 =(r 0 +s 0 )/2 

A partir da solupào aproximada obtida x 0 , fazemos a seguirne verificapào: 
Se f(x 0 ) = 0, entào x 0 è urna raiz de f(x) 

Se f{r 0 )f(x 0 ) < 0, entào r, = r 0 e s 1 = x 0 
Se f(r 0 )f(x 0 ) > 0, entào ^ = x 0 e s 1 = s 0 

Neste momento, temos o novo intervalo \r p s ; l, cuja amplitude é igual à 
metade da amplitude do intervalo [r 0 , s 0 ] e que contém a raiz desejada x. 

O procedimento é repetido novamente, isto é, calcula-se o ponto mèdio do 
intervalo [r p s r ] dado por x 1 = (r 7 +s 7 ) /2, conforme observado na Figura 3.8. 



Figura 3.8 Intervalo contendo o ponto mèdio. 

A partir da solupào aproximada obtida x p fazemos a seguirne verificapào: 
Se f(x 1 ) = 0, entào x 1 è urna raiz de f(x) 

Se f(r 1 )f(x 1 )<0, entào r 2 =r 1 e s 2 =x 1 
Se f(r 1 )f(x 1 )>0, entào r 2 =x 1 e s 2 =s 1 

Assim sucessivamente, temos um intervalo genèrico [a s ( ], em que calcu- 
la-se o ponto mèdio deste intervalo, dado por x j = (/;. + s, ) /2\ 

Se f(Xj) = 0, entào x t è urna raiz de f(x) 

Se f{r j )f(x j )< 0, entào r l+1 =r, e s i+1 =x i 

Se f(r) i f(x j )> 0, entào r M =x ; e s i+1 =s i 










Convergència 

Podemos observar que no Mètodo da Bissegào determinamos urna raiz 
da equagào, construindo sequèncias de intervalos g e s, e urna sequència de 
solugòes aproximadas x„/ = 0,1,... 

Essa sequència de solugòes aproximadas è convergente para a solugào 
desejada, urna vez que os intervalos sào divididos pelos pontos médios corres- 
pondentes e sào renomeados de forma que a raiz permanega dentro do interva- 
lo. Ver mais detalhes na Referència [1]. 

Estimativa do nùmero de iteragòes 

O nùmero de iteragòes necessàrias para se obter urna raiz x da equagào 
f(x) = 0, pelo Mètodo da Bissegào com urna precisào e > 0, previamente fixada, 
decorre do seguinte: 

Supondo-se que x està entre x n e s n , temos: 


X n- X \^( X n 



(s 0 -r 0 ) 

£n+1 


Impondo que 



para garantir que \x n -x 


<e temos: 


log 



\ 

< log (£) 


ou ainda, 


n> log(s 0 -r 0 )-logje) ? 
log 2 

Logo, n è o nùmero minimo de iteragòes que devem ser realizadas para 
obter x com urna precisào e. 

Algoritmo 3.1 

1. Dados e>0, o intervalo inicial [ r 0 ,s 0 ] que contenha a raiz, isto è, 
f{r 0 )f(s 0 ) < 0. Faga Pare = Falso, / = 0 

2. Enquanto Pare = Falso faga: 

2.1 Determine x, = (g +s),/2 











2.2 Se | f(Xj) | < £, entào Pare = Verdade 
Senio 

Se f{r i )f{x j )< 0, entào r i+1 =r) e s M =x, 
Senào r ì+1 = x, e s i+1 =s i 

\x. — X- 

2.3 Se — ^—r-^ < £ < e entào Pare = Verdade 

M 

Senào i = i + 1 


Exemplo 3.3 

Considere a seguirne fungào f{x) = 4x- e x . 

Usando o Mètodo da Bissegào, podemos determinar x tal que f(x) = 0 
com urna precisào e = 0.01, conforme segue: 

a) Inicialmente determinamos graficamente urna vizinhanga para a raiz, con- 
sideramos a forma equivalente 4x = e x , ou seja, f 1 (x) = 4x e f 2 (x) = e x , 
conforme ilustrado anteriormente na Figura 3.5. 

Observando a Figura 3.5, podemos concluir que a raiz x encontra-se na 
intersecgào dos gràficos /j(x) e f 2 (x), e pertence ao intervalo [0, 1\ 

b) Considerando o intervalo inicial r 0 =0 e s 0 =1, temos f{0)f{1)<0, pois 
f(0) = -1 ef(1) = 1.2817 e, portanto, temos realmente urna raiz no inter¬ 
valo [0,1]. Observe que a fungào f(x) = 4x-e x é continua no intervalo 
dado. 

c) Construindo a sequència de solugòes aproximadas, temos que 
x 0 = ( r 0 +s 0 )/2 = (0 + 1)/2 = 0.5000 é urna solugào inicial. 

Como f(r 0 )f(x 0 ) = (-1)(0.3513)<0, temos que o novo intervalo r ; e s 1 
serà dado por: 

l r i= r o = ° 

}s, = x 0 = 0.5 


Calculamos a nova solugào aproximada, conforme segue: 


x 1 = (r, +s 1 )/2 = (0 + 0.5 ) /2 = 0.2500 





Novamente, verificamos o critèrio de parada: 


x i~ x o 


= 1 > £, e corno este 


nào està satisfeito, repetimos o processo para calcular as novas solugòes apro- 
ximadas, corno segue: 

Como f(r 1 )f(x 1 ) = (-1)(-0.2840)>0, temos que o novo intervalo r 2 e s 2 
sera dado por: 


\r 2 =x 1 =0.25 
}s 2 = s r =0.5 


Calculamos o ponto mèdio do intervalo [r 2 , s 2 ] e temos: 
x 2 = (r 2 +s 2 )/2 = ( 0.25 + 0.5 ) /2 = 0.3750 


Verificamos novamente o critèrio de parada, dado por 



= 0.3333 > e, 


e corno este nào està satisfeito, sucessivamente repetimos o processo de càlculo 


das novas solugòes aproximadas, corno exibidas a seguir: 


x 3 = 0.3125 -» I X . 2 ' = 0.2000 > e 


x. = 0.3438 -» I * 3 ' = 0.0910 > e 

x 4 


x 5 = 0.3594 -> \ Xs * 4 I = 0.0434 > e 


x 6 = 0.3516 I Xe * 5 ' = 0.0222 > e 


x 7 = 0.3555 I X f Xg I = 0.0110 > £ 


x T 


xo = 0.3574 I X f X? I = 0.0053 < £ 

6 y 


Como o critèrio de parada està satisfeito, temos que a solugào aproximada 
de f(X) = 0 è dada por x = x 8 = 0.3574. 
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Podemos calcular o nùmero minimo de iteragòes usando a expressào: 


log(s 0 - r 0 ) -logje) 
log 2 


1 = 5.6439 


Portanto, n > 5.6439, isto é, devemos executar no minimo 6 iteragòes para 
obter a raiz x com a precisào e desejada, o que pode ser comprovado no exem- 
plo dado anteriormente com n = 7 iteragòes. 

3.3.2 Mètodo de Newton 

Seja f:3ì^3ì urna fungào continua e diferenciàvel, desejamos determinar 
a raiz ou raizes de f(x), ou seja, determinar x tal que f(x) = 0. 

O Mètodo de Newton consiste em, a partir de urna solugào inicial aproxima- 
da x 0 dada inicialmente, determinar urna sequència {x,}, / = 1,2,... de solugòes 
aproximadas para a raiz de f(x) da seguinte forma: 

Tragamos a reta tangente ao gràfico da fungào f(x) no ponto x ( , e no pon¬ 
to em que essa tangente cortar o eixo das abscissas, temos a primeira solugào 
aproximada x i+1 para a raiz x. 

Novamente no ponto x M tragamos a reta tangente ao gràfico de f(x), e no 
ponto em que essa nova tangente cortar o eixo das abscissas, temos a segunda 
solugào aproximada para a raiz, e assim sucessivamente até determinarmos, 
com urna tolerància prefixada e, a raiz x desejada. 

Ilustramos graficamente, conforme Figura 3.9: 
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Figura 3.9 Mètodo de Newton. 













Definindo corno a o àngulo formado com o eixo das abscissas através da 
reta tangente ao gràfico da fungào f(x) no ponto x, (Figura 3.9), temos: 

\_ f(*i) 

J (*/-*/+») 

Do Càlculo Diferencial Integrai, sabemos que tg(a) è a derivada da fungào 
f(x) no ponto x ( , isto é, tg(a) = f\x i ) e, assim, podemos escrever: 



f\x i ) = 


f(Xj) 

(Xi~X i+1 ) 


Portanto, temos o processo iterativo chamado de Mètodo de Newton, corno 
segue: 


x i+ i=x i - 


f(Xj) 

f'(Xi) 


i = 0,1,2,... 


Observagào 

O Mètodo de Newton è também conhecido corno Mètodo das tangentes, 
devido a sua interpretagào gràfica naturai. 

Convergència do Mètodo de Newton 

A convergència do Mètodo de Newton pode ser tratada usando o Teorema 
de convergència do Mètodo das Aproximagòes Sucessivas, o qual pode ser vis¬ 
to na Referència [1]. 

Definigào 3.2 Convergència quadràtica 

Dizemos que um mètodo iterativo possui convergència quadràtica se 

e ì+ i 

lim—^- = k, em que k è chamada constante assintótica de proporcionali- 

e f 

dade, e j =|x ( .-x| e e j+1 = | x i+1 -x| sào os erros cometidos nas iteragòes 
correspondentes. 

Teorema 3.1 

O Mètodo de Newton possui convergència quadràtica. 

Prova: Referència [1]. 

Para observar as boas propriedades de convergència do Mètodo de 
Newton, exibimos o seguinte exemplo: 
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Exemplo 3.4 

O càlculo do nùmero 43, usando o Mètodo de Newton, consiste na resolu- 
gào da equagào x 2 - 3 = 0 e, usando urna tolerància tixa com e = 10 4 , temos: 


*M = X i 


^1 = X, 

f (*/) 


(xf-3) 

(2x,) 


A partir de urna solugào x 0 inicial, geramos a sequència de solugòes 
aproximadas: 

x 0 = 1.5000 -» solugào inicial dada 
x, = 1.7500 -> I X \ ~ X ° I = 0.1429 > e 

\ x i I 

x P = 1.7321 -» \ \ = 0.0103 > e 

| X 2 | 

X q = 1.7321 ^l = 0 . 0QQ0<e 

xJ 


Como o critèrio de parada està satisfeito, podemos parar e observar que 
essa sequència converge para x = 1.7321 = 43. 

Dessa forma, podemos observar que à medida que os valores de x k apro- 
ximam da raiz x, a convergència torna-se muito ràpida, isso devido à proprieda- 
de da convergència quadràtica do Mètodo de Newton. 

Algoritmo 3.2 


1. Defina as fungòes f(x), f\x) e considere e > 0 urna tolerància fixa. 

2. Escolha x 0 urna solugào inicial. 

Faga Pare = Falso e / = 0 

3. Enquanto Pare = Falso faga: 

f(Xi) 


3.1 x j+1 = x, - 


f (*,) 


\ X — X \ 

3.2 Se J— 4 —entào Pare = Verdade 
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Senào i = i + 1 











Observagào 

O leitor pode incluir nesse algoritmo urna modificagào no critèrio de para¬ 
tìa, considerando o valor da fungào no ponto x h isto é, | f(x ( .) | < e. 

Exemplo 3.5 

Usando o Mètodo de Newton, resolvemos a equagào cos(x) - x = 0 com 
£ = 0 . 001 . 


A partir do processo iterativo x i+1 = x, 
x 0 = 0.7 solugàoinicialdada 
x. = 0.7394 I X ! X ° I = 0.0533 > e 

\ x i I 

x ? = 0.7394 \ \ = 0.0000 > e 

\x 0 


f(x j ) „ . 

, ' , geramos a sequencia: 

f(x,) 


Como o critèrio de paratìa està satisfeito, tomamos corno solugào aproxi- 
mada para f{x) a solugào x = x 2 = 0.7394. 

Observagào 

Podemos, ainda, modificar o Mètodo de Newton, da seguinte forma: 

O valor calculado da derivada na 1 a iteragào, f\x 0 ) = k, em que ke % è 
fixado e substituido no processo iterativo de Newton durante as iteragòes: 

Assim, temos: 


o qual è conhecido corno Mètodo Modificado de Newton. 

3.3.3 Mètodo das Secantes 


O Mètodo das Secantes consiste em modificarmos o Mètodo de Newton, 


cujo processo iterativo è dado por x j+1 = x, - 
fungào f (x,-) da seguinte forma: 


f(Xj) 

f'W' 


aproximando a derivada da 









f\x,) = 


(Xl-Xl-l) 


Observe que nesse caso, estamos trocando a inclinagào da reta tangente 
pela inclinagào da reta secante à curva, conforme pode ser observado Figura 
3.10 a seguir. 

Substituindo a expressào aproximada da derivada f\xj) no Mètodo de 
Newton, podemos escrever: 

f(Xj) 

Y — Y — ___ !— _ 

1+1 1 fjxj-fjx^) 

(*/”*/-») 


Simplificando a expressào anterior obtida, temos o Mètodo das Secantes 
conforme segue: 

X ' +1 f(x i )-f(x j _ 1 ) 

Assim, dados os pontos x i _ 1 e x ( , tragamos a reta secante passando por 
(Xj_ v f{x j _ 1 )) e {x j ,f{x i )) \ onde està cortar o eixo das abscissas temos a aproxi- 
magào x i+1 para a raiz x, conforme ilustrado na Figura 3.10. 



Figura 3.10 Mètodo das Secantes. 

Assim, na Figura 3.10, considerando o àngulo a formado pela reta secante 
à curva no eixo das abscissas, podemos escrever: 
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tga= rw = m-nx M ) 

(X,-X M ) I X - X, ,, 


Assim, temos o Mètodo das Secantes dado por: 

_ Xj-i f(Xj ) — X, f(Xj ? ) 

Algoritmo 3.3 

I.Seja f(x) continua e e>0 urna tolerànciafixa. 

2. Escolha x 0 e x 1 duas solugòes aproximadas iniciais. 
Faga Pare = Falso e / = 0 

3. Enquanto Pare = Falso fapa: 

31 x f(x l )-x j f(x l _ 1 ) 

1 + 1 

\ x. — x ■ \ 

3.2 Se J—< s, entào Pare = Verdade 

I X M I 

Senào i = i + 1 


Convergència 

Como o Mètodo das Secantes è urna modificapào do Mètodo de Newton, 
as condigòes de convergència sào parecidas, observando-se que nào temos 
mais a propriedade de convergència quadràtica. 

Quando f(x k ) = , podemos ter problemas de convergència, isto è, a 

sequència gerada pelo mètodo pode divergir. 

Exemplo 3.6 

Considere a equagào x 2 + 4x + 3 = 0] usando o Mètodo das Secantes com 
e = 0.01, podemos resolver essa equagào conforme segue: 

Podemos escrever a equagào dada na forma equivalente x 2 =-4x-3. 

Chamando f 1 (x) = x 2 e f 2 (x) = - 4x -3, temos urna vizinhanga para as 
raizes, na intersecgào dos gràficos de /j(x) e 400, conforme Figura 3.11 : 










Figura 3.11 Localizagào de rafzes. 


Usando o processo iterativo do Mètodo das Secantes 
y x jf( x i~i) e tornando corno solucòes iniciais x n =-0.5 e 

x 1 =-2.0, obtemos a solugào aproximada para a raiz, corno segue: 


Xgfjx^-x^jXg) _ (-0.5)(-1)-(-2)(1.25) _ 3 

f( Xl )-f(x 0 ) (-1) - {1.25) -2.25 


Calculamos o critèrio de parada, dado por: 



0.5000 > £ 


84 


Como o critèrio de parada nào està satisfeito, determinamos as demais 
solugòes aproximadas, corno seguem: 


x 3 = -0.5001 -» I X , * 2 ' = 1.6661 > £ 


x, = -1.0769 -» I = 0.5356>e 


x 5 = -1.0159 -» I X ? X . 4 ' = 0.0600 > £ 


x 6 = -0.9996 -» I ^ Xg I = 0.0163 > e 


Xy = -1.0000 -> I X | Xg I = 0.0004 < £ 
















Como o critèrio de parada està satisfeito, temos a solugào da equagào 
x = -1.0000. 

Observagào 

Podemos também determinar a outra raiz negativa, conforme exibida no grà¬ 
fico da Figura 3.10, apenas mudando as solugòes iniciais e novamente usando o 
processo iterativo do Mètodo das Secantes, até o critèrio de parada ser satisfeito. 

3.4 Exercfcios 

1. Usando os métodos da Bissegào, de Newton e das Secantes, resolva a 
equagào 3cos(x)-e-°- 4x =0 e e = 0.1. Localize graficamente urna vizinhanga 
para as raizes. 

2. Seja a fungào f(x) = e x ~ 2 +x 5 - 4. Determine o valor de x em que a fun- 
gào f(x) = 2, usando o Mètodo de Newton com e = 0.001. 

3. Determine x tal que f(x) = 0 , em que f(x) = x 3 - cos(x) + 2, usando: 

a) Mètodo da Bissegào com e = 0.1 

b) Mètodo de Newton com e= 0.001 

c) Mètodo das Secantes com e= 0.01. 

4. Determine um ponto da fungào f(x) = 6x 4 - 3x 3 -2x-3, em que a pri- 
meira derivada de f(x) se anula, usando o Mètodo de Newton com e =0.001. 
O ponto encontrado è um ponto de minimo da fungào dada? Justifique teorica¬ 
mente suas afirmagòes. 

5. Considere a fungào f(x) = 3x 2 + sen(x) - 2. Localize graficamente urna 
vizinhanga para as raizes de f(x). 

a) Determine a raiz negativa de f(x) usando o Mètodo de Newton com 
£ = 0 . 001 . 

b) Determine a raiz positiva usando o Mètodo das Secantes, com 
£ = 0 . 1 . 

6. Determine as raizes da fungào f{x) = 0.2x 3 -3.006x 2 + 15.06x - 25.15, 
usando um dos métodos vistos com e = 0.001. 

7. Considere a fungào f(x) = x-sen(x)-0.5. Determine x tal que f(x) = 0, 
usando: 

a) Mètodo das Secantes com e = 0.01. 




b) Mètodo de Bissegào corri e = 0.1. 

c) Determine urna raiz para f(x) dada, com o Software Numèrico, usando 
o Mètodo de Newton e e = 0.0001 - Referència [1]. 

8. A concentralo de bactérias poluentes C em um lago diminuì de acordo 
com a seguinte fungào: 


C(t) = 75e~ 15 ‘ + 20e~° ° 75t 


Deseja-se determinar o tempo necessàrio para que a concentragào de 
bactérias seja reduzida ao valor de 15. 

a) Resolva o problema usando o mètodo gràfico. 

b) Resolva o problema usando o Mètodo de Newton com aproximagào 
inicial t = 6 e e = 0.001. 

9. Usando o Software Matlab, piote o gràfico e determine urna raiz das fun- 
gòes abaixo: 

a) f(x) = 1 / [(x - 0.3) 2 + 0.01]+ [1/(x -0.9) 2 + 0.04]-6 

b) f(x) = e x -cos(x) -2 

c) f(x) = 2x 3 + ln(x) -5 

d) f(x) = x 5 -6.7x 4 + 8.4x 3 -10.8x 2 + 8x -6.8 

10. Faga um mapa conceitual detalhado sobre solugào numèrica de equa- 
gòes, considerando os métodos vistos, introduzindo labels indicando se aprendeu 
(A), nào aprendeu (N), se gostou (G), detestou (D ), achou interessante (I), etc. 




UNIDADE 4 


Interpolagào e Aproximagào de fungòes 



4.1 Introdugào 


Apresentamos nesta unidade a aproximagào de urna fungào de urna va- 
riàvel reai por outras fungòes mais simples, de modo que operagòes em geral 
sejam realizadas com mais facilidade. 

Essa situagào ocorre quando temos urna fungào f(x) e està apresenta 
um grau de dificuldade, por exemplo, para avaliar em pontos, derivar ou ainda 
integrar, ou mesmo quando conhecemos essa fungào em um nùmero finito de 
pontos de um intervalo [a, b], geralmente obtida em experimentos, sem o co- 
nhecimento de sua forma analitica. 

Dessa forma, essa fungào sera aproximada por fungòes polinomiais, expo- 
nenciais, trigonométricas, etc., que representarào a fungào originai, e para obter 
qualquer informagào sobre a fungào originai, utilizamos a sua forma aproximada. 

Exibimos nesta unidade apenas a aproximagào de urna fungào f(x), uti¬ 
lizando fungòes polinomiais, pela simplicidade no tratamento, continuidade e 
diferenciabilidade dessas fungòes, embora possamos estudar outros tipos de 
aproximagòes, corno as indicadas nas Referèncias. 

Posteriormente, apresentaremos o Mètodo dos Minimos Quadrados para 
aproximagào de urna fungào definida em um nùmero finito de pontos, e nesse caso, 
consideramos a minimizagào da soma dos quadrados dos erros nos pontos. 

Como ilustragào, apresentamos o problema da temperatura em um lago, 
conforme Figura 4.1. 



Figura 4.1 Lago. 

A temperatura T(°C) em um lago numa determinada regiào e època do ano 
è dada em fungào da profundidade P(m), conforme dados a seguir: 










T(°C) 

23.8 

22.8 

21.8 

20.8 

19.8 

18.7 

11.1 

P(m) 

0 

2.3 

4.9 

9.1 

13.7 

18.3 

22.9 


Deseja-se saber qual a profundidade do lago nas temperaturas 17.3°C e 
13.4°C, que nào foram medidas. 

Como tais valores nào se encontram na tabela, para resolver esse proble¬ 
ma, necessitamos do conhecimento de interpolagào, corno segue: 

4.2 Interpolagào Polinomial 

Considere umafungào f(x) definidaem x 0 ,x p ...,x n , n + 1 pontos distintos 
de um intervalo [a, b ], e denotamos y, = f(x j ), i = 0,1,..., n conforme represen- 
tagào na Figura 4.2: 



Figura 4.2 Valores da fungào f(x) nos pontos de interpolagào. 

Desejamos aproximar a fungào f(x) definidaem x 0 , x p ..., x n , n + 1 pontos 
distintos de um intervalo [a, b] por um polinòmio P(x) de grau menor ou igual 
n, tal que este coincida com a fungào nesses pontos, isto é, P(x i ) = f(x j ) = y,, 
i = 0,1,...,n. 

O Teorema seguinte garante a existència e a unicidade do polinòmio que 
desejamos determinar. 

Teorema 4.1 Existència e Unicidade 

Considere urna fungào f(x) definidaem x 0 ,x p ...,x n , n +1 pontos distintos 
de um intervalo [a, b], entào existe um ùnico polinòmio P(x) de grau menor ou 
igual a n, tal que P(x j ) = f(x i ) = y i , i = 1,...,n. 





















Prova: 


Seja um polinòmio de grau n, P(x) = a n x n + a n _ 1 x n 1 +.... + a 1 x + a 0 , tal 
que P(x,) = f(Xj) = y ; , / = 0,1,..., n. 

Dessa forma, para que o polinòmio coincida com a fungào nos n + 1 pon- 
tos, temos: 

Vo + a n -1 X V + 

a n X 1 + a n-1 X V + 

< 

ax n n + a n iX 1 }- 1 + ... + a.x + a n = y n 

L n n n-1 n 1 n 0 J n 


+ a,x 0 + a 0 = y 0 
+ a 1 x 1 +a 0 = y 1 


Podemos observar que obtemos um sistema de equapòes lineares Ax = b, 
em que x = (a n , a n _„..., aj, b = (y 0 ,y 1 ,...,yj e a matriz A dada por: 


*3 

v n-1 
*0 

.. X, 

x ? 

x ; n -' . 

.. X. 

X n 

v n-1 

.. X, 


O det(A), chamado de determinante de Vandermonde, é dado por 

cfef(/\) = n(x / -x y ) = (x,-x 2 )(x,-x 3 )...(x n _, -x n ) 

i<i 


Como os pontos x,,/' = 0,1,...,n sào distintos, podemos mostrar que 
det(A)#0, o que significa que o sistema linear possui urna ùnica solugào e, 
portanto, os coeficientes a 0 ,a p ...,a n do polinòmio sào unicos, dados pela reso- 
lugào desse sistema. Assim, o polinòmio P(x) existe e é ùnico. 

Definigào 4.1 Polinomio interpolador 

Denominamos polinòmio interpolador de urna fungào f(x) definida em 
x 0 ,x p ...,x n , n+1 pontos distintos de um intervalo [a, b] o polinòmio P(x) de 
grau menor ou igual a n, tal que coincide com a fungào nos pontos x (J / = n 
istoé, P(x j ) = f(x j ) = y j , i = 1,...,n. 

Representamos graficamente, conforme Figura 4.3: 








Figura 4.3 Interpolalo. 

Embora o polinòmio interpolador P(x) coincida com a fungào nos pontos 
de interpolalo x 0 , x p ..., x n , espera-se que P(x) = f(x) para x*x,, i = 0,...,n, 
ou seja, estimamos f(x) pelo polinòmio interpolador e cometemos um erro nessa 
aproximagào, dado por: 

E(x) = f(x)-P(x) 

Podemos representar, graficamente conforme Figura 4.4: 



Figura 4.4 Erro de interpolagào. 

Apresentamos a seguir urna expressào geral para o erro cometido quando 
aproximamos urna fungào f(x) por um polinòmio interpolador P(x). 






















Teorema 4.2 


Seja f(x) urna fungào definida em x 0 ,x p ...,x n , n + 1 pontos distintos de 
um intervalo [a, b] e n + 1 vezes diferenciàvel. Se P(x) interpola f(x) nesses 
pontos, entào o erro cometido em um ponto x é dado por: 

E(x) = f(x)-P(x) = -^-f^ n+1 \E > ) 

(n + 1)! 

n 

emque \\i(x) = Y\(x-x i ) e £e [. x 0 ,x n ]. 

i=0 

Prova: Referència [1]. 

4.3 Limitante superior para o erro 

Na expressào do erro do Teorema 4.2, o paràmetro £, nào é conhecido no in¬ 
tervalo [x 0 , x„] e, portanto, nào é possivel calcular o valor numèrico de f {n+1) (^). 

Dessa forma, apresentamos urna estimativa para o erro corno segue: 

E(x) = f(x)-P(x) = -^-f^ n+1 \E > ) 

(n + 1)! 

Podemos escrever: 


E(x) 


¥(*) 

(n + 1)! 


f in+1) (%) 


| vW | 

(n + 1)! 


f {n+1) (t,)\< l ¥(X) l M 
1 (n + 1)! 


com M= max||T n+? >(x)|, xe [x 0 ,x n ]|. 


Assim, temos um limitante superior para o erro: 


E(x) I < I ¥(X) I M. 


(n + 1)! 












Observagòes 

Podemos calcular urna estimativa para o erro somente quando tivermos a 
expressào analitica da fungào f(x), pois de acordo com a fòrmula do limitante 
superior para o erro, temos que dispor da (n +1) -ésima derivada dessa fungào. 

Nos casos em que tivermos apenas a fungào tabelada em um nùmero finito 
de pontos, sabemos que estamos cometendo um erro no ponto a ser avaliado, 
mas nào é possivel estimà-lo. 

Embora a resolugào do sistema linear obtido na prova do Teorema 4.1 for- 
nega urna maneira para determinar o polinòmio interpolador de urna fungào, apre- 
sentamos, também nesta unidade, a Fòrmula Interpolatória de Newton-Gregory. 
Nas Referèncias podemos ver alguns estudos sobre outras fórmulas interpolató- 
rias, corno a Fòrmula de Lagrange e a Fòrmula de Newton, as quais sào fórmu¬ 
las diferentes, mas que levam no mesmo polinòmio interpolador, pois corno foi 
demonstrado, este é ùnico. 

Para o desenvolvimento do polinòmio interpolador de Newton-Gregory, 
apresentamos o conceito de diferengas finitas, corno segue: 

4.4 Diferengas finitas 


Considereumafungào f(x) continua no intervalo [a, b], Sejam x 0 , x p ...,x n , 
n +1 pontos distintos deste intervalo [a, b ], tais que x j+1 - x l ■ = h, i = 0,1,... n-1, 
isto é, os pontos sào equidistantes. 

Definigào 4.2 Diferenga finita de ordem zero 

A diferenga finita de ordem zero de urna fungào f(x) definida nos pontos 
x e [a, b] é dada por: 

A °f(x) = f(x) 

Definigào 4.3 Diferenga finita de ordem n 

A diferenga dividida de ordem n de urna fungào f(x) definida nos pontos 
xe [a, b] é dada por: 


A n f(x) = A n 1 f(x + h) - A n 1 f(x) 

Nessa definigào, para r = 1 temos o operador diferenga finita progressivo 
A dado por: 
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A f(x) = f(x + h)-f(x) 


Assim, desenvolvendo-se os operadores podemos escrever 

A °f(x) = f(x) 

A 1 f(x) = f(x + h)-f(x) 

A 2 f(x) = f(x + 2h) - 2f(x + h) + f(x) 

A 3 f(x) = f(x + 3h) - 3f(x + 2h ) + 3f(x + h)- f(x) 







f n ì 


A n f(x) = 

f) 

f(x + nh)~ 


f(x + (n-1)h) + ... + (-1) n 


f(x) 


) 





Exemplo 4.1 

Considere urna fungào tabelada nos pontos, corno segue: 


x < 

0.1 

0.2 

0.3 

f(Xi ) 

6 

12 

15 


Càlculo das diferengas finitas de ordem zero: 

A °f(x 0 ) = f(x 0 ) = 6 
A °f(x 1 ) = f(x 1 ) = 12 
A °f{x 2 ) = f{x 2 ) = 15 

Càlculo das diferengas finitas de ordem 1 : 

a7(x 0 ) = A °f( Xl )- A °f(x 0 ) = (12-6) = 6 
A7(x ? ) = A °f(x 2 ) - A °f(x 1 ) = (15 -12) = 3 















Càlculo das diferengas finitas de ordem 2: 


A 2 f(x 0 ) = AV(x ; ) - A7(x 0 ) = (3 - 6) = - 3 

Podemos organizar o càlculo das diferengas finitas, conforme a Tabela 4.1 
a seguir: 

Tabela 4.1 Diferengas finitas. 

A°f A'f À 2 f À 3 f 


A°f(x 0 ) 

ATfXg) 


A°f(x i ) 


A 2 f(x 0 ) 


ATf^) 

A 3 f(x 0 ) 

A°f(x 2 ) 

AT(x 2 ) 

A 2 f (x i ) 

A°f(x 3 ) 




Os valores das diferengas finitas estào dispostos na tabela da seguirne forma: 

aV(x,) = A°f(x 1 ) - A°f(x 0 ) 

A 2 f(x 0 ) = A 1 f(x 1 )-A 1 f(x 0 ) 

A 3 f(x 0 ) = A 2 f(x 1 )-A 2 f(x 0 ) 

Exemplo 4.2 

Considere urna fungào f(x) tabelada nos pontos, corno segue: 



0.2 

0.4 

0.6 

0.8 

f(*i) 

7.8 

5.9 

10.8 

13.2 


Conforme Tabela 4.1, construfmos as diferengas finitas: 



M— 

CO 

<1 

H— 

OJ 

<1 

< 

o 

<1 

0.2 

7.8 


-1.90 

0.4 

5.9 6.80 


4.90 -9.30 

0.6 

10.8 -2.50 


2.40 

0.8 

13.2 
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4.5 Fórmulas Interpolatórias 

4.5.1 Fòrmula Interpolatola de Newton-Gregory 


Considere umafungào f(x) definida em um intervalo [a, b] e, x 0 , x p ..., x n , 
n + 1 pontos distintos e equidistantes deste intervalo. A Fòrmula interpoladora 
de Newton-Gregory é dada pelo seguirne Teorema 4.3: 

Teorema 4.3 

Seja urna fungào f(x) definida em x 0 ,x p ...,x n , n + 1 pontos distintos e 
equidistantes de um intervalo [a, b]. O polinòmio 

_ 0 , A ’f(x n ) A 2 f(x n ) 

P( x ) = A°f(x 0 ) + (x-x 0 ) + (x-x 0 )(x-x 1 ) 2jh ° + - + 

+ (X - x 0 )(X - X, ) ... (X - X n _-j ) A ^ X h ° n ) 


interpola a fungào f{x) dada nestes pontos, isto é, P(x,) = / = 0,1,... 

Prova: Referència [1]. 

Como desenvolvido anteriormente, urna estimativa para o erro é dada pelo 
limitante superior: 




| V(x) | 

(n + 1)! 


max 


f (n+1 \x) 


,xe [x 0 , x n ] 


Exemplo 4.3 

Considere a fungào f(x) = ln(x) tabelada nos pontos, corno segue: 


X, 

2 

3 

4 

f(X,) 

0.6931 

1.0986 

1.3863 


Podemos determinar o polinòmio interpolador Newton-Gregory, avaliar 
f(3.2) e calcular um limitante superior para o erro, corno segue: 

Usamos os très pontos tabelados para obter um polinòmio interpolador de 
grau <2: 


P(x) = A°f(x 0 ) + (x-x 0 ) 


A V(x 0 ) 
Uh 1 


+ (x 


x 0 )(x 


xù 


A 2 f(x 0 ) 

2 !h 2 
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Conforme Tabela 4.1, construfmos as diferengas finitas: 



H— 

OJ 

<1 

< 

H— 

O 

<1 

2 

0.6931 


0.4055 

3 

1.0986 -0.1178 


0.2877 

4 

1.3863 


Assim, temos o polinòmio interpolador de Newton-Gregory: 


P(x) = 0.6931 + (x-2) 0-4055 + (x - 2)(x - 3) 0/1178 


=-0.0589x 2 + 0.7000X -0.4713 


Portanto, f(3.2) = P(3.2)=1.1656 
Limitante Superior para o erro: 

I vi/(x) I , 

E(x) <-!- '-max A n+? >(x) xe[x n ,x 1 

1 (n + 1)! 1 1 0 n 


Como a fungào f(x) = ln(x), temos que as suas derivadas: 

fO){x) = —,x*0, f (2 \x) = x * 0 e A 3 >(x) = —. 
x x 2 x 3 


2 

Podemos observar que a fungào f (3 \x) = — é decrescente em mòdulo no 

x J 

intervalo [2,4], portanto, max | f 3 (x) | = 0.2500, em x = 2. em x = 2. 

Assim, 


E(3.2)< 


(■ 3.2-2){3.2-3){3.2-4) 
3! 


(0.25) = 0.0080. 


4.5.2 Polinòmio Interpolador de Newton-Gregory com mudanga de variàvel 


Podemos, ainda, escrever o polinòmio interpolador de Newton-Gregory, 
considerando a seguirne mudanga de variàvel: 












u = — —— ou x = x n +uh 
h 0 

Como os pontos sào equidistantes, x r =x 0 + rh segue que: 

(■ x-x r ) = (u-r)h 

Note que nessa nova variàvel u os pontos de interpolalo x 0 ,x 1 ,...,x n 
correspondem aos pontos 

u = 0,1,2,..., n. 

Dessa forma, podemos escrever o polinòmio de Newton-Gregory na variàvel u: 

A 1 f(x n ) A 2 f(x n ) 

P{u) = A °f(x 0 ) + (u-0) — + ( u ~ °^ u - 2! + ' ' ' 

A n f(x n ) 

+(u - 0){u — Y 1 

n! 


O limitante superior para o erro é dado por: 


E(u)< 


(n + 1) ! 


h n+1 M 


em que M = màx\ f (n+1 \x) , xe [a,b]. 

Exemplo 4.4 

Considere a fungào f{x ) da tabela nos pontos, conforme segue: 


*/ 

0.1 

0.2 

0.3 

f(Xi) 

2 

8 

10 


Podemos determinar o polinòmio interpolador de Newton-Gregory na varià¬ 
vel u e avaliar f(0.12) corno segue: 

O ponto x = 0.12 corresponde na variàvel u: 


u = 


(x-x 0 ) 


( 0 . 12 - 0 . 1 ) C2 


h 


0.1 
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. AV(x 0 ) _ A 2 f(x n ) 

P(u) = A°f(x 0 ) + (u -0) —^ + {u - 0)(u -1) —AyAZ 

Conforme Tabela 4.1, construimos a tabela de diferengas finitas: 



A°f A'f A 2 f 

0.1 

2 


6 

0.2 

8 -4 


2 

0.3 

10 


Assim temos: 


P{u) = 2 + u{6) + u{u-1)-^- = -2u 2 +8u + 2 


Portanto, f{0.12) = P u {0.2) = 3.52. 


4.6 Aproximagào de fungòes: Mètodo dos Minimos Quadrados 

O Mètodo dos Minimos Quadrados consiste em, dados os valores de urna 
fungào f(x) nos pontos x„ i = m , determinar urna fungào g(x) que melhor 
se aproxime dessa fungào conhecida. Para essa aproximagào, usaremos fun¬ 
gòes polinomiais, exponenciais, logaritmicas, trigonométricas, etc. 

Assim, conhecidos os valores de f(x) nos pontos x jt i = 1,..., m, desejamos 
determinar g(x), tal que g(x j ) = f(x i ), i = 

Como ilustragào, apresentamos o problema de previsào de escoamento de 
àgua, corno segue, conforme Figura 4.5: 









Figura 4.5 Escoamento de àgua. 

A estimativa dos tamanhos dos reservatórios depende de estimativas acu- 
radas do escoamento da àgua no rio que està sendo confinado. 

Para alguns rios, registros históricos de longa duragào de tais dados de 
escoamento sào dificeis de obter. Entretanto, os dados meteorológicos sobre 
precipitagào geralmente estào disponiveis mesmo após muitos anos. 

Portanto, é util determinar a relagào existente entre o escoamento e a pre¬ 
cipitagào. Essa relagào pode ser usada para fazer urna estimativa do escoamen¬ 
to nos anos em que apenas as medidas de precipitagào foram feitas. 

Os dados disponiveis para um rio que deve ser estancado estào dispostos 
na seguinte tabela: 


Precipitagào (cm) 

88.9 

108.5 

104.1 

139.7 

127 

94 

116.8 

99.1 

Escoamento (m 3 /s) 

14.6 

16.7 

15.3 

23.2 

19.5 

16.1 

18.1 

16.6 


Deseja-se determinar urna previsào do escoamento anual de àgua se a 
precipitagào for de 120 cm. 

Caso discreto 

Consideremos o caso em que a fungào f(x) é definida em um conjunto 
discreto, isto é, a fungào é conhecida em m pontos, geralmente obtidos em expe- 
rimentos, conforme a tabela: 


x i 


x 2 



f(x,) 

f(Xì) 

f(x 2 ) 


f(x m ) 





















Graficamente, temos a disposilo dos pontos obtidos no experimento, con¬ 
forme Figura 4.6: 



Figura 4.6 Dados de um experimento. 

Observando a disposigào dos pontos (x /; f(x,)), / = 1,..., m na Figura 4.6, 
vemos que g(x) possui o comportamento de urna reta, isto é, um polinòmio de 
grau 1 : 


g(x) = a 7 g 7 (x) + a 2 g 2 (x) = a,x + a 2 


Com g 7 (x) = x e g 2 (x) = 1. Assim, escolhemos urna famflia de fungòes as 
quais dependem dos paràmetros a 1 e a 2 . 

O problema agora consiste em determinar os paràmetros a 1 e a 2 de modo 
que a fungào g{x) melhor se ajuste aos dados da tabela. 

Para falar em melhor ajuste, temos que ter um critèrio para a escolha dos 
paràmetros a 1 e a 2 , isto é, temos que ter urna medida para o erro cometido 
nessa aproximagào. 

Definicào 4.4 Erro ou desvio 

Definimos e(x, ) = f(x,-) - g(x, ) corno o erro ou desvio, cometido numa apro- 
ximagào de urna funpào f{x) por urna fungào g(x), nos pontos x ; , / = m. 

Dessa forma, desejamos determinar urna fungào g(x) de modo que nos pon¬ 
tos x„ i = 1,..., m os desvios sejam pequenos. Nesse caso é tentador desejar que 

m 

a soma dos erros sejam mfnimos, isto é, que X e ( x /) se J a Entretanto, 

i=i 
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esse fato nào traduz que g(x) seja urna boa aproximagào para a fungào f(x), pois 

m 

podemos ter (x, ) = 0 , serri que os erros sejam pequenos. 

i=i 

m 

Poderiamos, também, considerar ]T|e(x ( .)| m i nima ’ porém esse critèrio 

i=i 

acarreta dificuldades de resolugào, pois a fungào valor absoluto nào é diferencià- 
vel na origem. 

Urna maneira para contornar esses problemas consiste em considerar urna 
medida para o erro da seguirne forma: 

m m 

Minimizar ^ e(x,) 2 = Minimizar (f(x,■) - g(x : )) 2 

i=1 hi 

Assim, considerando o exemplo da Figura 4.6, desejamos encontrar urna 
fungào g(x) = a 1 x + a 2 que melhor se aproxime da fungào f(x), de forma que 

m 

E(a p a 2 ) = ^(e(x,)) 2 seja minimo. Do Càlculo Diferencial de fungòes de vàrias 

hi 

variàveis, se a fungào E(a p a 2 ) possui um ponto de minimo, entào suas deriva- 
das parciais devem ser nulas nesse ponto, isto é: 


dE 

da 1 


= 0 


dE n 
e - = 0 

da , 


Derivando E(a v a 2 ) com relagào à variàvel a 1 , temos: 


dE_ = _d_ 
da 1 da 1 


ì=i 


da 1 


£(e(x,.)) 2 =— £(a r x,. +a 2 -f{x j )f 


hi 


= X 2{a 1 x j + a 2 -f(x j ))x j 


hi 


= 2 


a,£xf+a 2 £x ; 


hi 


hi 


-2 


IVI',) 


hi 


= 0 


Assim, temos: 


f m 2 'N 


( m 

X*, 

a 1 + 

X x , 



l hi J 


hi 


i) 

















Derivando E(a p a 2 ) com relagào à variàvel a 2 , temos: 


dE 

da 2 


c) 

da 2 



(e(x,)) 2 


\ 


v 


d 

da 2 


r m 




+ 3.2 


f(x,)) 2 

y 


= X 2{a lXj + 


i=1 


a 2 




= 2 


m m 

a /X x , +m(a 2 )-^f(x i ) 


i=1 


i=1 


= 0 


Assim, temos; 

m 

a 1 +ma 2 = Y J f{x i ) 

ì=i 

Portanto, os paràmetros a, e a 2 que minimizam o erro E(a p a 2 ) necessa¬ 
riamente satisfazem o seguirne sistema de equagòes lineares: 



< 


f m 2 N 


' m ^ 

X*, 

a, + 

l X i 

l ,=? J 


U - * J 

r m ^ 



l x i 

a 1 + ma 2 = 


a 2 = X 


X f ( x i) 

i=1 


X, f(x,) 


O sistema de equagòes obtido é chamado sistema de equagòes normais, o 
qual pode ser resolvido por qualquer mètodo numèrico visto anteriormente. 

Exemplo 4.5 

Considere a seguinte fungào f(x) da tabela corno segue: 


X i 

0.5 

0.8 

0.9 

1.0 

1.1 

f(x,) 

1.51 

1.82 

1.91 

2.00 

2.10 


Piote os pontos dados e verifique se os mesmos possuem um comporta¬ 
mento linear, isto è, de urna reta (verifique isso!). 

Usando o mètodo dos Minimos Quadrados, determinamos dentre todas as 
retas g(x) = a r x + a 2 aquela que melhor se ajusta aos dados. 





























5 5 

I*, 2 X*, 

i=1 i=1 


>' 


X^ x /) x / 

i=1 

tx> 5 

i=i 


_ d 2 _ 


-1 

_ 1 


Temos a seguinte tabela de dados: 



x, 

f(Xi ) 

x? 

f(x i )x j 

E 

0.5 

1.51 

0.25 

0.7550 

0.8 

1.82 

0.64 

1.4560 

0.9 

1.91 

0.81 

1.7190 

1.0 

2.00 

1.00 

2.0000 

1.1 

2.10 

1.21 

2.3100 

4.30 

9.34 

3.91 

8.2400 


Assim, temos o seguinte sistema de equagòes lineares: 


'3.91 

4.30 ' 

a i 


'8.24' 

4.30 

5 

_ 3 2 _ 


9.34_ 


Usando o Mètodo de Eliminagào de Gauss, na forma pivotada, temos: 

3.91 4.30 . 8.24 

0 0.2711 . 0.2785_ 

Assim, 

a 2 =1.0272 
a 1 = 0.9778 


Portanto, a fungào linear que melhor se aproxima dos dados da tabela é 
dada por g(x) = 0.9778x + 1.0272. 

Càlculo do erro: 

efx,) 2 = (f(0)-g(0)) 2 =0.0000 
e{x 2 ) 2 = {f(1)-g(1)) 2 =0.0001 
e(x 3 ) 2 = (f(2) - g(2)) 2 = 0.0000 
e(x 4 ) 2 = (f (3)-g(3)) 2 =0.0001 
e(x 5 ) 2 = (f (4) -g(4)) 2 =0.0000 
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5 

Portanto, ]T e(x,) 2 = 0.0001 e qualquer outra reta possui a soma dos qua- 

i=1 

drados dos erros superior a esse valor obtido. 

Podemos ter dados experimentais em que seja necessàrio aproximar a 
fungào f(x) por um polinòmio de grau 2, isto é, urna paràbola: 

g(x) = a l g 1 (x) + a 2 g 2 {x) + a 3 g 3 (x ) = a r x 2 +a 2 x + a 3 


com g 1 (x) = x 2 ,g 2 (x) = x eg 3 (x) = 1. 

Generalizando esse procedimento, escrevemos g(x) corno urna combina- 
gào linear de fungòes corno segue: 


g(x) = a,g r (x) + a 2 g 2 {x) +... + a n g n {x) 


com as fungòes g,(x) escolhidas. 

Procedendo de maneira anàloga ao caso do ajuste linear, a reta, podemos de- 

m 

terminar os paràmetros a j ,i = 1,...,n de forma que o erro ^e(x,) 2 seja minimo. 

i=i 

Assim, temos: 


dE < m 
— = 0 o 
da 1 


i=i 


a i + ^9 2 ( x i)9i( x i)) 


+ 


y v'= ? 

A m 


a 2 +. 


ì=i 


)0,(*, ) a n='L f ( X i)9 1 ( X i) 


i=1 


dE _ ( m 

= 0 <=> 


da 0 


i=i 


J,9 1 ( x ì)9 2 ( x ì) a i + ^9 2 ( x ì)9 2 ( x ì) 


+ 


y v ,w 

A m 


a 2 +. 


ì=i 


)0 2 (* ; ) a n='L f ( X i )9 2 ( X i) 


i=1 


d E n 1 m 
- = 0 o 


3a„ 


A f 




i=i 


i=i 


y v' 

f m \ m 

)0 n (* ( ) a n =Xf(x,.)g n (x,.) 


a 2 +... 


ì=i 


y 


ì=i 


+ 
























Portanto, para determinar os paràmetros a„i = 1 ...,n temos que resolver 
o seguirne sistema de equagòes lineares: 


A f 


^9 1 ( x l )9 1 ( x i) a,+ J\g 2 (x i )g 1 (x i ) 


y v 


a 2 +... 


f m 


'Lg n (x i )g 1 ( x i) a n = 'Lf(x i )g 1 (x i ) 

i=1 i=1 


f rn 


'Zg 1 ( x i)g 2 ( x i ) 


f m 




Z gJx, )g 2 ( x i ) a n = £ f ( x i)g 2 ( x i ) 


i=1 


i=1 


t,9l( X i)9n( X i) 
v ,=7 y 


a, + 


f m \ 

y Lg 2 ( x i)g n ( x i) 


a 2 +... 


f m 


i=1 i=1 


O sistema de equagòes lineares obtido é denominado sistema de equagòes 
normais, o qual pode ser resolvido por qualquer mètodo visto anteriormente. 

Resolvido esse sistema, determinamos os paràmetros a,, / = e, con- 
sequentemente, a fungào: 


g(x) = a 1 g 1 (x) + ... + a n g n (x) 


que melhor se ajusta à fungào f{x) nos pontos x p ...,x m no sentido dos mini- 
mos quadrados. 

Exemplo 4.6 

Considere urna fungào f(x) definida conforme tabela: 


*/ 

-2 

-1 

0 

1 

2 

3 

f( x i) 

19.01 

3.99 

-1.00 

4.01 

18.99 

45.00 


Dispondo os pontos num gràfico podemos ver que g{x) pode ser ajustada 
por urna paràbola. Verifique està afirmagào! 

Assim, tomamos g(x) = a 1 x 2 +a 2 x + a 3 , isto é, g 1 (x) = x 2 , g 2 (x) = x, 
g 3 (x) = 1 e determinamos os paràmetros a p a 2 e a 3 de modo que g(x) se ajus- 
te aos dados da tabela no senso dos minimos quadrados. 


107 































Temos o sistema de equagòes normais: 


6 6 6 

xv x *, 3 x *, 2 

i=1 i=1 i=1 

c c c 


V 


i=1 

ODO 

È *, 3 É x / 2 È*/ 

i=1 i=1 i=1 


a 2 

= 

i=1 

hi 2 t*i 6 

_ i=i i=i 


a 3_ 


S^) 

. i=i 


Temos a seguinte tabela de dados: 




x? 

*? 

xf 

f(Xi) 

x, f(Xi) 

xff(Xi) 


-2 

4 

-8 

16 

19.01 

- 38.02 

76.04 


-1 

1 

-1 

1 

3.99 

- 3.99 

3.99 


0 

0 

0 

0 

- 1.00 

0.00 

0.00 


1 

1 

1 

1 

4.01 

4.01 

4.01 


2 

4 

8 

16 

18.99 

37.98 

75.96 


3 

9 

27 

81 

45.00 

135.00 

405.00 

3 

19 

27 

115 

90.00 

134.98 

565.00 


Assim, podemos escrever o sistema de equagòes normais, corno segue: 


~115 

27 

19 

V 


565.00 

27 

19 

3 

a 2 

= 

134.98 

19 

3 

6 

a 3_ 


90.00 


Usando o Mètodo de Eliminagào de Gauss, na forma pivotada, temos: 

~115 27 19 . 565.00 

0 12.6604 -1.4612 . 134.98 

0 0 2.6926 . 90.00 

Assim: 

a 1 = 5.0893 
< a 2 = 0.0515 
a 3 =-1.1403 
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Portanto, g(x) = 5.0893x 2 + 0.0515 x-1.1403 


































Càlculo do erro: 


e(x r ) 2 = (f(x 1 ) -giXj)) 2 =0.0108 
e(x 2 ) 2 = (f(x 2 ) -g(x 2 )) 2 = 0.0086 
e(x 3 ) 2 =(f(x 3 ) -g(x 3 )) 2 =0.0197 
e(x 4 ) 2 = (f(x 4 ) -g(x 4 )) 2 = 0.0001 
e(x 5 ) 2 = (f(x 5 ) -g(x 5 )) 2 = 0.1088 
e(x 6 ) 2 = ( f(x 6 ) -g(x 6 )) 2 = 0.0332 

6 

Portanto, e(x,) 2 =0.1812 e qualquer outra paràbola ajustada possuem 

i=i 

a soma dos quadrados dos erros superiores a esse valor obtido. 

4.7 Exercfcios 

1. Seja urna fungào f(x) definida conforme a tabela abaixo: 


X 

0 

1 

2 

3 

f(x) 

-5.6 

-4 

8.1 

9 


Determine o polinòmio interpolador de f(x), avalie f(2.3) usando a Fòrmula 
Interpolatola de Newton-Gregory. 

2. Seja urna fungào f(x) definida conforme a tabela abaixo: 


X 

0.1 

0.2 

0.3 

0.4 

f(x) 

8 

3.9 

4.8 

-8.9 


Determine o polinòmio interpolador de f(x) nos pontos dados, usando a 
Fòrmula Interpolatola de Newton-Gregory na variàvel u. Avalie f(0.35). 

3. Estude o polinòmio interpolador de Lagrange e o limitante superior para 
o erro - Referència [1]. 

4. Seja a fungào f(x) = e x + cos(x) +1 tabelada corno segue: 


X 

0 

0.5 

1.1 

f(x) 

3 

3.53 

4.46 


Usando o polinòmio interpolador de Lagrange, avalie f(0,6) e um limitante 


superior para o erro. 
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5.0 calor especffico da àgua em fungào da temperatura em °C é dado por: 


t 

30 

35 

40 

45 

c 

0.99826 

0.99818 

0.99828 

0.99890 


Com o Software Numèrico - Referència [1] -, usando o polinòmio interpola- 
dor de Newton-Gregory, calcule aproximadamente o calor especffico para 

t = 37.5 °C. 

6. A partir dos seus conhecimentos de aproximagào de fungòes, resolva o 
problema da temperatura em relagào à profundidade de um lago, proposto no 
infoio deste capftulo, usando um polinòmio interpolador de grau 2. Analise os 
resultados obtidos. 

7. Usando o Mètodo dos Mfnimos Quadrados, determine g(x) = a 1 x + a 2 
que melhor se ajusta aos dados da tabela abaixo: 


X 

-2 

-1 

0 

1 

2 

3 

4 

f(x) 

-4.98 

-3.00 

-1.01 

0.99 

3.01 

4.98 

7.01 


Calcule ^e(x ( ) 2 e analise os resultados obtidos. 

i=1 

8. Com o Software Numèrico - Referència [1 ] -, e usando o Mètodo dos Mf- 
nimos Quadrados, determine urna fungào g(x) que melhor se ajusta aos dados 
da tabela abaixo: 


X 

-2 

-1 

0 

0.1 

1 

2 

3 

f(x) 

6.00 

1.01 

0 

0.12 

3.00 

9.98 

20.98 


7 

Calcule ^,e(x ( ) 2 e analise os resultados obtidos. 

i=i 

9. Considera o problema de previsào de escoamento de àgua, exibido neste 
capftulo: 

a) Faga um gràfico dos dados da tabela. 

b) Usando o Mètodo dos Mfnimos Quadrados, determine a melhor reta 
que melhor se ajusta aos dados da tabela e faga urna previsào do es¬ 
coamento anual de àgua se a precipitagào for de 120 cm. 

10. Usando o Software Matlab, resolva os exercfcios 7) e 8). 




























11 .Faga um Mapa Conceitual detalhado sobre os tópicos de interpolagào, 
introduzindo labels indicando se aprendeu (A), nào aprendeu (N), se gostou (G), 
detestou (D), achou interessante (I), etc. 

12. Faga um Mapa Conceitual detalhado sobre os tópicos de ajuste de cur- 
vas - MMQ -, introduzindo labels indicando se aprendeu (A), nào aprendeu (N), 
se gostou (G), detestou (D), achou interessante (I), etc. 





UNIDADE 5 


Integragào Numèrica 



5.1 Introdugào 


Nesta unidade, apresentamos alguns métodos numéricos para calcular 
aproximadamente o valor da integrai de urna fungào com urna variàvel reai defi- 
nida num intervalo fechado [a, b]. 

De maneira geral, termos: 

b 

1 = J f(x)dx 

a 

Em que a fungào f(x) é continua com derivadas continuas no intervalo 
[a, b]. 

Sabemos que o valor da integrai I é dado pelo Teorema Fundamental do 
Càlculo para integrais corno: 

b 

I = J f(x)dx = F(b) - F(a) 

a 

em que F(x), chamada de fungào primitiva de f(x) é tal que F\x) = f(x). 

b 

Graficamente, podemos representar I = J f(x)dx corno a àrea A, formada 

a 

entre a curva e o eixo das abscissas quando f(x) assume valores positivos em 
[a, b], conforme a Figura 5.1 a seguir: 



Figura 5.1 Àrea sob o gràfico da fungào f(x). 

Exemplo 5.1 

' x 6 , 

Para calcular \x dx, corno F(x) = — satisfaz F(x) = x , entào: 
i 6 
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J x 5 dx = F(1) - F(0 ) = —- = 1/6 

o 6 6 

Muitas vezes, o càlculo da fungaci primitiva torna-se algo trabalhoso, em al- 
guns casos conhecemos apenas os valores da fungào tabelados em um nùmero 
finito de pontos e nào podemos calcular a integrai dessa fungào usando apenas 
os conhecimentos do Càlculo Diferencial Integrai, corno citamos anteriormente. 

Métodos numéricos sào desenvolvidos para calcular aproximadamente o 

b 

valor da integrai I = J f{x)dx, em que é necessàrio apenas o conhecimento da 

a 

fungào f{x) em um nùmero finito de pontos. Esse processo é conhecido corno 

Fórmulas de Quadratura de Newton-Cotes. 

Para ilustrar o tòpico de integragào numèrica, consideremos o problema de 
estimar medidas das margens até um rio conforme Figura 5.2: 



A partir de urna linha reta próxima às margens de um rio, foram feitas me¬ 
didas em (m) entre essa linha reta e as margens do rio, de 15 em 15 metros, a 
partir do ponto tornado corno origem. Tais dados foram registrados conforme a 
tabela abaixo. 



0 

15 

30 

45 

60 

M, 

0 

1 

0 

2 

4 

m 2 

5 

7 

7 

8 

10 


Deseja-se determinar o valor aproximado da àrea coberta pelo rio. 















5.2 Integragào numèrica usando interpolagli) 


Este mètodo consiste em aproximar urna fungào f{x) definida em 
x 0 ,x p ...,x n , n + 1 pontos distintos do intervalo [a, b], pelo polinòmio interpola- 
dor P n (x), e integrar o polinòmio P n (x) ao invés da fungào f(x). 

Assim, temos: 

X X 

J f (x)c/x = J P n (x)dx 


Podemos representar graficamente, conforme Figura 5.3: 



Figura 5.3 Integragào numèrica. 


5.3 Fórmulas de Quadratura de Newton-Cotes 


Considere urna fungào definida em x 0 , x p ...,x n , n + 1 pontos distintos e 
equidistantes do intervalo [a,b], isto é, x (+; -x, =/7, i = 0,...,n-1, com h>0 a 
distància entre os pontos. 

Nesse caso, podemos considerar o polinòmio interpoladorde Newton-Gregory 
para f(x)\ 


P(x)=D°f(x 0 ) + (x-x 0 ) 


DVjx,,) 

Uh 


+ (x-x 0 )(x-x 1 ) 


D 2 f(x 0 ) 
2 ih 2 


+ ... + (x-x 0 )...(x-x n _ ? ) 


A n f(x n ) 


n!h n 













Dessa forma, temos: 


jf(x)dx= J P n (x)dx 


Como f(x)=P n (x) + E n (x) em que o erro é dado por: 


( n+1) 

EJx) = (x-x„)...(x-x n )AAL-, x 0 <%<x 


Temos 


J f{x)dx = J P n (x)dx + J E n (x)dx 


Como os pontos sào distintos e equidistantes, podemos considerar o poli¬ 
nòmio de Newton-Gregory na variàvel u, em que: 

u = (x - x 0 )/h 
Assim, 


P„{u) = A°f(x 0 ) + uA 1 f(x 0 ) + U(U 2! — A 2 f(x 0 ) +. 

+ u(u-1)...(u-(n-1)) 

n! °' 


e o erro na interpolalo na variàvel u é dado por: 


= u(u-1)...(u-n) 

" (n+1)! 0 " 


Dessa forma, o erro na integrapào pode ser escrito corno: 


E n = \ E n {x)dx = h] E n (u)du = h\ 


u{u-1)...(u-n) 
(n + 1)! 


(n+1) 


h n+1 f(x) du 









Erro cometido na Integracào Numèrica 


X, 


X, 


Podemos aproximar J f(x)dx por j P n (x)dx e riessa aproximagào come- 

x o x o 

temos um erro expresso pelo seguirne teorema: 

Teorema 5.1 

a) Se f(x) possui n +1 derivadas contmuas no intervalo [ x 0 , x n ], e os pon- 
tos x, = x 0 + ih, i = 0,1,...,n subdividem o intervalo de integragào num 
nùmero fmpar de intervalos, entào existe um ponto t ,, tal que: 



0 


b)Se f(x) possui n + 2 derivadas continuas no intervalo [x 0 ,x n ], eospon- 
tos x, =x 0 + ih, i = 0,1,...,n subdividem o intervalo de integragào num 
nùmero par de intervalos, entào existe um ponto £, tal que: 



Prova: Referència [1]. 

5.3.1 Regra dos Trapézios 

Considere urna fungào f(x) definida em dois pontos x 0 e x 1 no intervalo 
[a, b]. O polinòmio interpolador da fungào f(x) de grau n = 1 é dado por: 



e 


X, 


J f{x)dx = J P 1 (x)dx = h jP,(u)du 


X, 


0 


em que u = ——^ 
H h 
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Representamos graficamente, conforme Figura 5.4: 



Figura 5.4 Regra dos Trapézios. 


Assim, 


"7 ' 

J f(x)dx = hj (A °f{x 0 ) + utif(x 0 )) du 


= hj A °f{x 0 )du + h j i/A 7 f(x 0 )du 


Como A °f(x 0 ) = f(x 0 ) e A 1 f{x 0 ) = f{x 0 )-f(x 0 ), temos 


7 ' ' 

J f(x)dx = /ij A°f{x 0 )du + hj uA 1 f(x 0 )du 


= hf(x 0 )u 


h (fiXj ) — f(x Q )) 


u 




hf(x 0 )+-(f(x 1 )-f(x 0 )) 


:-(f(x 0 ) + f(x,)) 



























Portanto, termos: 


j f(x)dx = ^[f(x 0 ) + f(x 1 )] 


denominada Regra dos Trapézios. 

Observagào 

O càlculo da integrai J f{X)dx corresponde à àrea do trapèzio formada 

Xo 

entre o polinòmio interpolador P ? (x) e o eixo das abscissas, conforme a Figura 
5.4, justificando, assim, a denominagào Regra dos Trapézios. 

Erro na Regra dos Trapézios 

Nesse caso, corno o intervalo de integragào n = 1 é impar, e usamos a parte 
(a) do Teorema 5.1, temos: 


—j(u)(u-1)du, x 0 <^<x 1 


i 

Como ju(u-1)du 

0 


—, temos que o erro é dado por: 




12 


f (2) & 


Limitante superior para o erro 

Podemos observar que o argumento £, na fòrmula do erro, nào é urna 
grandeza numèrica conhecida no intervalo [x 0 , x y ] e, portanto, nào è possivel 
calcular o erro cometido na Regra dos Trapézios. Dessa forma, podemos calcu- 
lar urna estimativa para o erro, que serà o limitante superior para o erro: 


Ei 


12 




12 


P 2 >(^e[x 0 ,x ; ]. 


Como | P 2) (^) | < maxU T 2) (x)|,xe [x 0 ,x r ]j-, temos que: 


E 1 |<^màx{| P 2 >(x) |,x 0 <x<x ? } 
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Exemplo 5.2 


2 

Podemos calcular J (sen(x) + x)dx usando a Regra dos Trapézios e um limi- 

i 

tante superior para o erro da seguinte forma: 

Tabelando a fungào f(x) nos pontos x 0 =1.0 e x 1 = 2, temos: 


X 

1 

2 

f(x,) 

1.8415 

2.9093 


Usando a Regra dos Trapézios, temos: 

2 h 1 

J (sen(x) + x)dx = — [f(x 0 ) + f(x ? )] = —[1.8415 + 2.9093] = 2.3754 


Limitante superior para o erro: 

E 1 \ < — màx{\ A 2, (x)|> 7 < x < 2} 


( 2 ) 

Como afungào f(x) = -sen(x), entào 


f(2) 



= max {| T 2 >(x) \,1<x<2} = 1 


pois a fungào | f (2 \x) | = sen(x) assume o valor màximo em x = 

Dessa forma, temos: 

i i (I ) 3 

E , <T^(7) = 0.0833 
in 

5.3.2 Regra dos Trapézios Generalizada 

A Regra dos Trapézios Generalizada consiste em dividirmos o intervalo 
de integragào [a, b] em n subintervalos iguais, de amplitude h, de forma que 
x 0 =a e x n =b e, aplicarmos a Regra dos Trapézios em cada um deles, isto é, 
em cada dois pontos consecutivos. 
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Graficamente, temos: 



Figura 5.5 Regra dos Trapézios Generalizada. 

x " h h 

J f(x)dx = -[f(x 0 ) + f(x 1 )] + - [f(x 1 ) + f(x 2 )] + ... 


+^lf(x n - 1 ) + f(x n )] = 

= — [f(x 0 ) + 2f(x 1 ) + 2f(x s ) + ... + 2f(x n _ 1 ) + f(x n )] 


Erro total na Regra dos Trapézios Generalizada 

O erro total na fòrmula generalizada é obtido a partir da soma dos erros 
cometidos em cada subintervalo, isto é: 

Como vimos em cada aplicagào da Regra dos Trapézios, temos a seguinte 
expressào para o erro: 


Assim, o erro total cometido é dado por: 

n tj 3 (2) 
i=1 1 ^ 

Como x j _ 1 < <Xj, i = 1 , ri e f {2 \x ) é urna fungào continua por hipóte- 
se, entào existe £ e [x 0 ,x n ], tal que 
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X f (2) (^/ ) = nf (2) (£)* x 0 -£>- x n 

i=1 

Assim, a expressào para o erro na Regra dos Trapézios Generalizada 
torna-se: 

h 3 (2) 

E T=-jj nf (&^[ x o’ x nl 

Como o nùmero de subintervalos n é dado por n = (x n -x 0 )/h, temos: 

h 2 (2) 

E t =-—(x n - x 0 )f(^),^e[x 0 , x n ] 

Limitante superior para o erro 

Como o argumento £ nào é conhecido, nào podemos determinar o erro 
exatamente, mas podemos calcular um limitante superior para o erro. 

t?2 

| E t \^J^( X n~ X o) màx {\ f (2) (X)\,X 0 <X<X n } 


Exemplo 5.3 

Usando a Regra dos Trapézios Generalizada, podemos calcular o valor 

i 

aproximado da integrai J(x/n(7 + x))dx com 5 subintervalos e um limitante su¬ 
ri 

perior para o erro. 

Solugào: 


n = 



—> 5 = 


( 1 - 0 ) 

h 


->h = 0.2 


Dessa forma, podemos tabelar a fungào f(x) corno segue: 


*1 

0 

0.2 

0.4 

0.6 

0.8 

1.0 

f(Xi) 

0 

0.0365 

0.1346 

0.2820 

0.4702 

0.6931 














Assim, temos: 


7 f_ 

J (x liti + x))dx= — [f(x 0 ) + 2(f(x 1 ) + f(x 2 ) + f(x 3 ) + f(x 4 )) + f(x 5 )] = 

= ^-(2.5397) = 0.2538 

Limitante Superior para o erro: 

| E t \^^ X n~ X o) màX {\ f (2) PO|,* 0 <X<X n } 

(2 + xì 

Como a funcào f (2) =--/ é urna funcào decrescente em mòdulo em 

(7+x) 2 

[0 /] (verifique isso, observando que a fungào f (3 \x ) < 0 no intervalo de integra¬ 
fo), entào màx| A 2 >(x) | = 2 emx = 0. 

Assim, 

E t < ^(x n - x 0 ) màx{| A 2 >(x) |, 0 < x < f}= ^(7 - 0)(2) = 0.0333 

5.3.3 Regra 1/3 de Simpson 

Considere urna fungào f(x) definida em très pontos distintos x 0 , x 1 e x 2 
e equidistantes no intervalo [a, b] e o polinòmio interpolador da fungào f(x), de 
grau n = 2, com x 0 = a e x 2 = b. 

Assim, 

_ A 7(x„) A 2 f(x n ) 

P 2 (x) = A°f(x 0 ) + (x-x 0 ) 0 +(x-x 0 )(x-x r ) ° 

e 

*2 *2 2 

J f(x)dx = j P 2 {x)dx = hjP 2 {u)du 

x o x o 0 


em que u = 


(*-* 0 ) 


h 
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Graficamente, temos: 



Figura 5.6 Regra 1/3 de Simpson. 
Assim, 


J f{x)dx = hj f(A °f{x 0 ) + uA 1 f(x 0 ) + U(U 1 7) A 2 f{x 0 ) 

x 0 0 V ^ * 

2 2 

= hj A °f{x 0 )du + /?J uA 1 f(x 0 )du 

O 0 

+ h 2 j^_ll & 2 f(Xii)d u 
0 * ■ 

Sabendo-se que 

A f(x 0 ) = f(x 0 ), A J f(x 0 ) = f(x 1 ) - f(x 0 ) 

A 2 f{x 0 ) = f(x 2 ) + 2f(x 1 ) + f(x 0 ) 

Temos: 

x 2 2 2 2 

J f(x)dx = h j A 0 f(x 0 )du + h J uAU{x 0 )du + h J 


du = 




= hf{x 0 ){u) 


/7(/(X ; )-/(X 0 )) 


f u ^ 


k 2 j 


0 






















L. 

- [f(x 0 ) + 4f(x 1 ) + f(x 2 )] 

U 


Assim, obtemos: 


J f(x)dx = ^[f(x 0 ) + 4f(x 1 ) + f(x 2 )] 


denominada Regra 1/3 de Simpson. 

Erro na Regra 1/3 de Simpson 

Nesse caso, o intervalo de integragào foi subdividido em um nùmero par 
de subintervalos, n = 2k, k = 1, 2,3,... e portanto, pela parte (b) do Teorema 5.1, 
temos: 

( n+2) 



0 


Assim, para n = 2 temos: 



o 



0 



Como o argumento £ nào é conhecido, nào é possivel calcular o erro exa- 
tamente, portanto, trabalhamos com um limitante superior para o erro: 
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Limitante Superior para o erro 

Temos que o erro é dado por: 


e 2 = - 


( 4 ) 

h 5 f& 

90 


<x<x„ 


Tornando em mòdulo a expressào dada anteriormente, temos que: 


E 2 = 


h*_ 

90 


f i4) &) Ue [x 0 ,x 2 \ 


Dessa forma, podemos escrever o limitante superior para o erro corno: 


h 5 


E 2 \ <—{màx\ n 4 \x) \,x 0 <x<x 2 } 


Exemplo 5.4 

2 

Usando a Regra 1/3 de Simpson, podemos calcular je x dx e um limitante 
superior para o erro da seguirne forma: 0 

Tabelando a fungào f(x) com h = 1.0, temos: 


x i 

0 

1.0 

2.0 

f( x i) 

1 

2.7183 

7.3891 


Dessa forma, temos: 

2 

J f(x) dx = - [f(x 0 ) + 4f{x 1 ) + f(x 2 )] 

0 J 

= -[1 + 4(2.7183 ) + 7.3891] = 6.4208 
3 

Portanto, 


2 

J e x dx = 6.4208 

o 












Limitante superior para o erro 

Da fòrmula do limitante superior para erro, temos que: 



h5 

<— {maxi A 4 >(x)I 
90 1 1 1 


xe[0,2 ]} 


Como a fungào f (4) (x) = e x , entào 


| f {4) (2.0) | = max{| A 4 >(x )\,xs [0,2]}= 7.3891, 

pois a fungào | A 4) (x) | é crescente no intervalo [ 0, 2\ 
Assim, 

I E. \<^7.3891 = 0.0821 
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5.3.4 Regra 1/3 de Simpson Generalizada 

A Regra 1/3 de Simpson Generalizada consiste em dividirmos o intervalo 
[a, b ] de integragào em n subintervalos iguais de amplitude h, em que n é um 
nùmero par de subintervalos, de forma que x 0 = a, x n = b , e aplicarmos a Regra 
1/3 de Simpson em cada 2 subintervalos, isto é, a cada 3 pontos consecutivos. 

Graficamente, temos: 



Assim, aplicando a Regra 1/3 de Simpson para cada 3 pontos, isto é, a 
cada 2 subintervalos, temos: 















j f(x)dx = ~ [f(x 0 ) + 4f(x 1 ) + f(x 2 )] + ^ \f(x 2 ) + 4f(x 3 ) + f(x 4 )] + ... 

*o J J 

+^[f(x n . 2 ) + 4f(x n _i) + f(x n )] 

= | [f(x 0 ) + 4(f(x 1 ) + f{x 3 ) + ... + f(x n1 )) 

+2{f{x 2 ) + f(x 4 ) + ...+ f(x n _ 2 ) + f(x n )] 

Erro total na Regra 1/3 de Simpson Generalizada 

O erro total cometido na Regra 1/3 de Simpson Generalizada é obtido a 
partir da soma dos erros cometidos a cada aplicagào da Regra 1/3 de Simpson 
em dois subintervalos subsequentes, isto é: 

Para cada dois subintervalos, 3 pontos distintos e equidistantes, aplicamos 
a Regra 1/3 de Simpson e o erro cometido para cada aplicagào é dado por: 

( 4 ) 

h 5 f(L) 

E 2 = -* 2,-2 * *>, z x 2 i’ i = 1,2,..., n/2 

Assim, temos: 

n/2 U5 (4) 

f( ^ì )> x 2i-s * li ^ x 2i , i = 1,2,..., n/2 

ì=i yu 

Como x 0 < '%! < x n , i = 1, 2 ,..., n/2 e a fungào f {4 \x ) é continua por hipóte- 
se, entào existe um [x 0 , x n \ , tal que: 

n/2 (4) 

)= n fm,x 0 <^< x n 


n/2 (4) 


Assim, a expressào para o Erro na Regra 1/3 de Simpson Generalizada é 
dada por: 


h 5 n <4) 

E t=-gÒ2 mX o^ X n 






(x — X ) 

Como o nùmero de subintervalos n = — n 0 , temos que: 


M ( 4 ) 

E f =—(x n -x 0 )f(^),x 0 <^<x n 
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Limitante Superior para o erro 

Como [x 0 ,x n ] nào é urna grandeza numèrica conhecida, trabalhamos 
com um limitante para o erro, conforme segue: 

h 4 

\ Et \-J80 {Xn ~ X 0 >)màX {\ f(4)(X) l’ X 0- X - X n} 


Exemplo 5.5 

Usando a Regra 1/3 de Simpson Generalizada, 

f 1 

J- dx com 4 subintervalos e um limitante superior para 

Q (l + X) 

forma: 


podemos 
o erro, da 


calcular 

seguinte 


Temos h = (Xf1 Xp) =- = 0.5 
n 4 


Tabelando a fungào f(x) com h = 0.5, temos: 


X, 

0 

0.5 

1.0 

1.5 

2.0 

f(Xi) 

1 

0.6666 

0.5000 

0.4000 

0.3333 


Temos 

2 j #_ 

J -- [f(x 0 ) + 4f(x 1 ) + 2/(x 2 ) + 4f(x 3 ) + f(x 4 )] = 

0 (7 + X) o 

= —[7 + 4((0.6666 ) + 0.4000) 

3 

+2(0.5000) + 0.33333] = 1.1000 


Assim, 


J 


_ 1 _ _ 

(1 + x) 


dx = 1.1000 


















Limitante Superior para o erro 

Temos 


I ~W) {Xn ~ Xo) màx i\ fi4)(x) I’ x °- x - x "} 


180 


24 

Como f { \x) = -- é urna funcào decrescente em mòdulo no intervalo 

( 1 + x ) 5 

[0,2],màx\f( 4 \x)\ = 24, em x=0. 

Assim, 

0 5 4 

\E t < ^-{2 - 0)(24) = 0.0166 
! M 18Q 

5.4 Exercicios 


0.7 

1. Calcule I = J (e 6 * + 8.4x + 3)dx usando: 

o.i 

a) Regra dos Trapézios com 6 pontos. 

b) Regra 1/3 de Simpson com 6 subintervalos. 

c) Calcule um limitante superior para o erro em cada caso. 


2. Seja urna fungào f(x) tabelada corno segue: 


X 

-1 

0 

1 

2 

3 

4 

5 

f(x) 

1 

0.51 

0.42 

0.82 

1.91 

0.99 

1.88 


5 

Calcule J f(x)dx usando a Regra dos Trapézios e a Regra 1/3 de Simpson. 

i 

3. Calcule j(4.6x 2 +cos(x)) + 3.5dx usando a Regra dos Trapézios com 6 

o 

pontos e um limitante superior para o erro. 

4. Estude a Regra 3/8 de Simpson e um limitante superior para o erro e 
resolva o exercicio 2) - Referència [1]. 

2 

5. Calcular J(e 3 * +8.5x + 6)dx usando: 

i 


a) Regra dos Trapézios com 4 subintervalos. 

















b) Regra 1/3 de Simpson com 6 subintervalos. 

c) Regra 3/8 de Simpson com 9 subintervalos e um limitante superior para 
o erro. 

6. Determine o menor nùmero de subintervalos em que podemos dividir o 
intervalo [0 1] para obter 

i 

I = j (5 e- 2x )dx 
0 

usando a Regra dos Trapézios com erro menor ou igual a 0.001. 

7. Resolva o exercicio 6) usando a Regra 1/3 de Simpson, a Regra 3/8 de 
Simpson e para a divisào encontrada calcular a integrai I. 


8. Urna linha reta foi tragada de modo a tangenciar as margens de um rio 
nos pontos A e B. Para medir a àrea entre o rio e a reta AB, foram tragadas per- 
pendiculares em relagào a AB com um intervalo de 0.05 m. Qual é essa àrea? 


Perpendiculares 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Comprimento(m) 

3.28 

4.02 

4.64 

5.26 

4.98 

3.62 

3.82 

3.68 

3.26 


Faga urna interpretagào gràfica do problema e resolva-o usando seus co- 
nhecimentos de Càlculo Numèrico. 

9. Com o Software Matlab, usando a Regra dos Trapézios, calcule as inte- 
grais das seguintes fungòes: 

5 

a) j(sen(x)-x 2 + 3.1)dx 

i 

3 

b) j (e x +cos(x)- 2.8)dx 

10. Com o Software Numèrico - Referència [1] -, usando a Regra 1/3 de 
Simpson, calcule o valor aproximado da integrai das seguintes fungòes: 

4 

a) J( e x + 5x + 4.9)dx , usando 6 subintervalos e um limitante superior 

i 

para o erro. 

2.5 

b) J (x 3 +6x + 2)dx , usando 9 pontos e um limitante superior para o 


erro. 















II.Faga um Mapa Conceitual detalhado sobre métodos numéricos para 
Integragào Numèrica, introduzindo labels indicando se aprendeu (A), nào apren- 
deu (N), se gostou (G), detestou (D ), achou interessante (I), etc. 




UNIDADE 6 


Solugào Numèrica de Equagòes Diferenciais 
Ordinàrias e Parciais 



6.1 Introdugào 


O estudo de Equagòes Diferenciais deve-se ao fato de que, em situagòes 
pràticas nas diversas àreas cientfficas, frequentemente temos problemas cuja 
modelagem matemàtica é representada por urna equagào diferencial. 

O desenvolvimento de métodos numéricos para resolugào dessas equa¬ 
gòes diferenciais é necessàrio, pois solugòes analiticas nem sempre sào possi- 
veis de serem determinadas. 

Tipicamente, problemas que envolvem urna fungào de urna variàvel e suas 
derivadas levam a urna Equagào Diferencial Ordinària, e problemas que envol¬ 
vem urna fungào de vàrias variàveis e suas derivadas parciais nos levam a urna 
equagào diferencial parcial. 

Um exemplo simples de equagào diferencial ordinària é o problema da pro- 
pagagào da fumaga de cigarro num ambiente fechado ou semifechado. 

Exemplo 6.1 

Considere urna sala de trabalho que contenha um volume V = 36 m 3 de ar 
originalmente isento de monóxido de carbono. Suponha que, a partir do instante 
t = 0, a fumaga do cigarro contendo 4% de monóxido de carbono é introduzida 
na sala por um fumante, com urna vazào de 0.002 m 3 por minuto, e que essa 
sala possui urna janela ou porta semiaberta, de forma que a mistura gasosa, 
distribuida homogeneamente, sai da sala com a mesma vazào. 

a) Qual é a expressào para a quantidade q(t) de monóxido de carbono na 
sala num instante t? 

b) Sabendo-se que a exposigào prolongada de monóxido de carbono com 
concentragào acima de 0.00012 é prejudicial à saude, em quanto tempo 
essa concentragào serà atingida? 


Para equacionarmos o problema, consideremos a quantidade q(t) de mo¬ 
nóxido de carbono na sala no instante t. A taxa de variagào de q(t) é igual a quan¬ 
tidade de monóxido de carbono que està entrando menos a que està saindo. 


dq 

dt 


- 0.002 

100 


0.002 

36 


<7(0 


A solugào analitica desse problema com a condigào inicial q(0) = 0 é a fungào 


200 200 - 

q(t) = -e 

139 139 


139 t 


2500000 
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Considerando = temos: 

36 


139 t 

f(t) = 0.04- 0.04e 2500000 

que é a porcentagem da quantidade de monóxido de carbono na sala no instante 
t, corno exibido na Figura 6.1 : 



Figura 6.1 Monóxido de carbono numa sala de fumantes. 

Para sabermos o instante em que a sala atinge o limite insuportàvel de 
monóxido de carbono, resolvemos a seguinte equagào: 

f(t) = & = 0.00012 
36 

E obtemos o instante t = 54.1 minutos, o que significa que em 54 minutos 
a sala estarà saturada de monóxido de carbono. 

Entretanto, muitas equagòes diferenciais nào possuem solugào analitica, o 
que nos leva ao desenvolvimento de métodos numéricos para obter a solugào 
aproximada. 

Definicào 6.1 

Urna Equagào Diferencial Ordinària de ordem n é urna equagào da seguinte 
forma: 


F(t,y(t),y'(t),y" (f y (n) (f)) =0 










em que estào envolvidas a fungào incògnita y = y(f) e suas derivadas até a 
ordem n, sendo t a variàvel independente. 

A notagào y {>) representa a derivada de ordem j da fungào incògnita y em 

dW 

relagào a variàvel independente t e pode também ser representada por 

Exemplo 6.1 

As equagòes a seguir podem ser escritas conforme a definigào anterior. Sào 
equagòes diferenciais ordinàrias, em que a fungào desconhecida é y = y(x)e 

dv d 2 v 

as suas derivadas sào y' = —,y" =——. 

7 dx dx 2 

dv 

a) — = 2x - y é urna equagào diferencial de ordem 1. 
dx 

b) y"+(1- y 2 )y '+ 6y = 0, em que y = y(t ) e suas derivadas y \t) = — 

dt 

d 2 v 

e y"(t) = , e urna equagào diferencial de ordem 2. 

6.2 Equagào diferencial de primeira ordem com valor inicial 

Um problema de valor inicial consiste em urna equagào diferencial 
y' = f(x,y), e urna condigào inicial y (x 0 ) = y 0 , em que y 0 é um valor da fungào 
y dado no ponto x 0 . 

Assim, temos um problema de valor inicial denotado por: 

\y' = f(x,y) 

0 ) 

{y(x 0 )=y 0 

Resolver o problema (1) consiste em determinar urna fungào y(x) que 
satisfaga (1), isto é, y' = f(x,y) e y(x 0 ) = y 0 . 

Exemplo 6.2 

Considere o problema do crescimento da poluigào numa cidade e o modelo 
matemàtico em que o crescimento da poluigào està sujeito à hipótese de que a 
taxa de variagào da poluigào p(t) é proporcional à poluigào num instante t, dado 
pela equagào diferencial: 

^ = kp , oup'(t) = kp(t) 


em que k é a taxa de variagào da poluigào anual. 
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Suponha que no instante t = 0 temos urna poluigào inicial de Ippm (urna 
parte por milhào). Se k = 0.065 por ano, podemos utilizar métodos analiticos ou 
numéricos adequados de resolugào e fazer a previsào dessa poluigào num perio¬ 
do de tempo t = 10 anos, resolvendo o problema de valor inicial: 


p' = 0.065 p 
P(0) = 1 


A solugào analitica dessa equagào diferencial é dada pela fungào 
p(t) = e 0,065t , pois ao substituirmos p(t) e sua derivada na equagào dada, veri- 
ficamos que a equagào està satisfeita. 

Urna condigào que garante a existència e a unicidade da solugào do proble¬ 
ma de valor inicial (1 ) é dada pelo seguirne Teorema 6.1 : 

Teorema 6.1 

Considere urna fungào reai f(x, y) continua no intervalo a<x<b, com a e 
b finitos e — o° < y < + oo. Se existe urna constante L tal que para todo xg [a, b] 
e para todo par de valores (y, y ; ) tivermos 

| f(x,y)-f(x,y 1 ) | <L\ y-y 1 \ 

conhecida corno Condigào de Lipschitz, entào existe urna ùnica fungào y = y(x) 
satisfazendo as seguintes condigòes: 

a) y(x) é continua e diferenciàvel para todo xg [a,b] ; 

b) y' = f(x,y(x)) para xs[a,b]\ 

c ) y(x 0 ) = y 0 , em que y 0 é um valor conhecido. 

Prova: Referència [11] 

Observagào 

Podemos observar pelo Teorema 6.1 que a solugào de urna equagào dife¬ 
rencial ordinària é urna fungào y = y(x) continua e diferenciàvel que satisfaz a 
equagào diferencial y' = f(x,y) e passa pelo ponto (x 0 ,y 0 ). 




6.3 Discretizagào 


O processo de discretizagào consiste em resolver numericamente urna 
equagào diferencial calciando aproximagòes para a fungào y = y(x) em pon- 
tos discretos x 0 ,x v x 2 ,...,x N do intervalo de discretizagào [a,b]. 

Para discretizar o intervalo [a, b], tomamos N subintervalos, N>1 , e faze- 

X — X 

mos x n =x 0 +nh, n = 0,1,2,..., N com x 0 =a e x N =b, sendo h= N N ° ■ 

Conhecendo y(x) num ponto inicial dado y(x 0 ) = y 0 (valor inicial), calcu- 
lamos passo a passo, nos pontos x 1 =x 0 + h, x 2 =x 0 +2h, x 3 =x 0 +3h,..., 
x n =x 0 +nh, as solugòes aproximadas y n para a solugào exata y(x n ), 
n = 0,1,2,...,N, conforme ilustrado na Figura 6.2: 



Figura 6.2 Discretizagào de um intervalo. 

O erro locai, cometido nas aproximagòes em cada ponto, é a diferenga en- 
tre o valor exato da equagào diferencial e o valor numèrico aproximado em cada 
um dos pontos do intervalo [a, b], isto é: 


e(x n ) = y(x n )-y n ,n = 1,...,N 



















6.4 Métodos Numéricos para Equagoes Diferenciais Ordinàrias 

6.4.1 Mètodo de Euler 

Seja a equagào diferencial de primeira ordem com valor inicial: 

ÌY' = f(x,y) 

i y(x 0 ) = y 0 

O Mètodo de Euler é baseado na expansào da fungào derivàvel y(x) em 
Sèrie de Taylor, isto è, expandindo a fungào y(x) por Sèrie de Taylor nas vizi- 
nhangas do ponto x n , até a ordem 1, podemos escrever: 

y(x n + h) = y(x n ) + hy\ x n ) 

Dessa forma, corno o ponto x n + h = x n+1 e y'(x n ) = f(x n ,y(x n )), temos: 
y(x n+1 ) = y(x n ) + hf(x n ,y(x n )) 


Assim, podemos escrever para n = 0: 


y(x,) = y(x 0 ) +W (x 0 ,y(x 0 )) 

e 


y(^) = y 0 + hf(x 0 ,y 0 ) 


Chamando y* =y(x 0 ) + hf(x 0 , y(x 0 )), temos o valor aproximado da fungào 
y(x) no ponto x p isto è, y(x 1 ) = y 1 . 

Portanto, o Mètodo de Euler no primeiro passo calcula y 1 =y 0 + hf(x 0 ,y 0 ) 
que è a aproximagào y 1 da solugào exata y(x ? ) no ponto x f . 

Repetindo esse procedimento para o ponto x 2 , temos no segundo passo 
y 2 =y 1 + hf {x p y r ) que è a aproximagào y 2 da solugào exata y(x 2 ) no ponto 
x 2 - E assim sucessivamente para cada um dos pontos x n , n = 0,1,2,...,N de 
um intervalo [a,b]. 




De maneira geral, temos: 


yn + 1=yn +h ^ X n’Vn)’ n = 0, l 2,... 

Podemos interpretar graficamente o Mètodo de Euler aproximando a fun- 
gào y(x) pela reta tangente ao gràfico de y no ponto dado ( x 0 , y(x 0 )), da se¬ 
guirne forma: 


(y-y 0 ) = m(x-x 0 ) 


em que m = y'(x 0 ) = f(x 0 ,y 0 ) é o coeficiente angular dessa reta. 
Assim, podemos escrever: 

( y-y 0 ) = f ( x o’yo)( x - x o ) y=y 0 +f( x 0 ,y 0 )( x - x 0 ) 


Tornando y = y 1 e x = x p temos que: 


yi =yo +f ( x o’y 0 )( x i- x o) 


Ou ainda 

yi=y 0 +hf( x o,y 0 ) 

Assim, o ponto (x 1 ,y 1 ) pertence a essa reta tangente, e o mesmo racioci- 
nio pode ser utilizado para os demais pontos discretizados, conforme ilustramos 
na Figura 6.3: 



Figura 6.3 Mètodo de Euler. 
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O erro locai cometido é a diferenga entre o valor exato e o valor aproximado 
em cada ponto (. x n ,y n ),n = 0,1,...,N. 


Erro de Troncamento 


por: E = 


No Mètodo de Euler, truncamos a Sèrie de Taylor em p = 1, e o erro è dado 

y"m 2 


2 ! 


[a b\. 


Como o paràmetro £ nào è conhecido no intervalo de discretizagào, vamos 
considerar urna estimativa para o erro, isto è, um limitante superior para o erro: 


E< 


Mh 2 

2 ! 


M = max {|y "(x)|, x e [a, £>]) 


Algoritmo 6.1 

Considere a equagào diferencial, com valor inicial 


J y' = f(x,y) 
\y(x 0 ) = y 0 


1. Declare 


a) Fungào f(x,y) 

b) Condigòes iniciais: y(x 0 ) = y 0 

c) Intervalo [a,b], em que a = x 0 

_ Q 

d) Nùmero de subintervalos N e calcule h = - 

' N 

2. Para n = 0,...,N -1, faga: 
infoio 

x n+1 =x n + h 

y n+ i=yn + hf ( x n’y n ) 

firn 


Exemplo 6.3 

Considere a seguinte equagào diferencial de ordem 1 com valor inicial: 


y' = f(x,y) = y 

y(x 0 ) = y(0) = 1 








a) Usando o Mètodo de Euler, podemos calcular a solugào aproximada 
y(1.0 ) com h = 0.2. 

b) Compare os valores exatos dessa equagào, com os valores aproxima- 
dos obtidos em a). 

Considerando h = 0.2, temos a discretizagào do intervalo dada por: 
x 0 = 0, x 1 = 0.2, x 2 = 0.4, x 3 = 0.6, x 4 = 0.8 ex 5 =1.0 

Do Mètodo de Euler: 

yn + 1=yn +hf (Xn,y„) 

Càlculo de y. 

Para n = 0, temos a condigào inicial y(x 0 ) = y(0) = 1, entào: 

Yi = y 0 + hf (x 0 ,y 0 ) 

= y 0 + 0.2 (y 0 ) = 1 +0.2(1) = 1.2000 

Càlculo de y 2 

Para n = 1, temos: 

y 2 = yi + hf(x 1 , y ; ) 

= y, + 0.2 (y 1 ) = 1.2000 + 0.2(1.2000 ) = 1.4400 

Càlculo de y 3 

Para n = 2, temos: 

y 4 = y 3 +hf(x 3 ,y 3 ) 

= y 3 + 0.2(y 3 ) = 1.7208 + 0.2(1.7208) = 2.0736 

Càlculo de y 4 


Para n = 3, temos: 




y 4 = y 3 +hf(x 3 ,y 3 ) 


= y 3 + 0.2(y 3 ) = 1.7208 + 0.2 ( 1.7208 ) = 2.0736 


Càlculo de y 5 

Para n = 4, termos: 


y 5 = y 4 + hf(x 4 ,y 4 ) 

= y 4 + 0.2(y 4 ) = 2.0736 + 0.2(2.0736 ) = 2.4883 


Portanto, termos que y(1) = y 1 = 2.4883. 

Estimativa para o erro 


E<—màx {| y"(x) |,xe [0,/]} 


Como a soluqào analitica y(x) da equagào é dada por y(x) = e x , temos que: 
max{| y"(x) |,xe [0,1]\= 2.7183, em x = 1, urna vez que y"(x) = e x é urna 
funqào crescente em mòdulo. 

Portanto, temos: 

E < { °' 2 ^ (2.7183) = 0.0544. 


Exibimos na Tabela 6.1 os valores exatos, a partir da soluqào exata 
y(x) = e x , a soluqào numèrica pelo Mètodo de Euler e os respectivos erros lo- 
cais nos pontos discretos. 


Tabela 6.1 


N 

*/ 

Sol. Exata 

y(*n) 

Sol. Aprox. 

y n 

Erro 

\y(x n )-y n \ 

0 

0 

1.0000 

1.0000 

0.0000 

1 

0.2 

1.2214 

1.2000 

0.0214 

2 

0.4 

1.4918 

1.4400 

0.0518 

3 

0.6 

1.8221 

1.7280 

0.0941 

4 

0.8 

2.2255 

2.0736 

0.1519 

5 

1.0 

2.7183 

2.4883 

0.2300 
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Observe na Tabela 6.1 que temos urna boa estimativa para o erro no infoio 
do processo, mas piora quando os pontos afastam-se do ponto inicial. 

Podemos observar também que para valores de h menores que h = 0.2 , 
obteremos melhores resultados, pois o erro depende da amplitude do passo h. 


6.4.2 Mètodo de Euler Aperfeigoado 


O Mètodo de Euler Aperfeigoado é urna modificagào do Mètodo de Euler 
desenvolvido anteriormente, da seguirne forma: 

1 ^ 

y n+ i =y n + ^ f ( x n’y n )+ f ( x n +^ /„ + h y n h n= 0 , 1 , 2 ,... 

Para melhor entendimento da expressào do Mètodo de Euler Aperfeigoado, 
faremos urna interpretagào gràfica, conforme Figura 6.4 a seguir: 



Figura 6.4 Mètodo de Euler aperfeigoado. 

Observando a Figura 6.4, vemos que: 

a) A reta t passa pelo ponto A = (x n ,y n ) e possui coeficiente angular 
y n = f(x n ,y n ). Usando o Mètodo de Euler, calcula-se y^ +1 =y n +hy' n \ 

b) A reta r passa pelo ponto B = (x n+1 ,y^ +1 ) e possui coeficiente angular 

f(x n +h,y n +hy' n ) = f(x n+v y E n+1 y, 

c) A reta r 1 passa pelo ponto B e sua inclinagào è a mèdia das inclinagòes 
das retas ter; 


d)A reta r 2 passa pelo ponto A e è paralela a reta ry 
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e) O valor y n+1 obtido pelo Mètodo de Euler Aperfeigoado é urna aproxima- 
gào para a solugào y(x) no ponto x n+1 . 


Algoritmo 6.2 

Considere a equagào diferencial 


\y\x) = f{x,y) 

ly(^ 0 ) = y 0 


1. Declare: 


a) Fungào f{x, y) 

b) Condigòes iniciais: y(x 0 ) = y 0 

c) Intervalo [a, b], em que a = x 0 

d) Numero de subintervalos N e calcule: h = - 

' N 

2.Para n = 0,1,2,3,...,N -1 

Calcule: 

infoio 

X n+ 1 = X n +h 
k 1 = f ( X n’Yn) 
k 2 =f( x n + i,y n + hk i) 

y n+ i =y n + ^ k i +k 2 ) 

firn 

Exemplo 6.4 

Considere a equagào diferencial com valor inicial: 


y' = f(x,y) = y 


y(x 0 ) = y(0) = 1 


Podemos, usando o Mètodo de Euler Aperfeigoado, determinar a solugào 
aproximada y{1.0) com h = 0.2 da seguinte forma: 

Considerando h = 0.2, temos a discretizagào do intervalo: 
x 0 = 0, x 1 = 0.2, x 2 = 0.4,x 3 = 0.6, x 4 = 0.8 ex 5 =1.0 






Càlculo de y, 

Para n = 0 temos a condigào inicial y(x 0 ) = y(0) = 1, entào 
k^fiXcYo) -+k 1 =y 0 = 1 

k 2 = f(x 0 + h,y 0 + hk 1 ) = y 0 + hk 1 =1 + 0.2(1) = 1.2000 

Do Mètodo de Euler Aperfeigoado temos: 

1 ^ 

yn + 1=yn+-^(k 1 + k 2 ) 

Portanto: 

y 1 = y 0 +^(k 1 +k 2 ) = 1 + ^~ (1 + 1.200) = 1.2200 

Càlculo de y 2 

Para n = 1 

k 1 =f(x 1 ,y 1 ) ->k 1 =y 1 =1.2200 

k 2 = f(x 1 + h,y 1 + hk 1 ) = (y 1 + hk 1 ) = 1.2200 + 0.2(1.2200) 

= 1.4640 

y 2 =y 1 +^(k 1 +k 2 ) = 1.2200+ ^-(1.2200+ 1.4640) = 1.4884 

Càlculo de y 3 

Para n = 2 

k 1 =f(x 2 ,y 2 ) -+k 1 =y 2 = 1.4884 

k 2 = f(x 2 +h,y 2 + hk 1 ) = (y 2 + hk 1 ) = 1.4884 + 0.2(1.4884) 

= 1.7861 

y 3 =y 2 +^(k 1 +k 2 ) = 1.4884 + (1.4884 + 1.7861) = 1.8159 

Càlculo de y 4 


Para n = 3 




k 1 = f ( x 3 ’Y 3 ) -» k 1 =y 3 = 1.8159 


k 2 =f( x 3 + h,y 3 + hk 1 ) = (y 3 + hk 1 ) = 1.8159 + 0.2(1.8159) 

= 2.1791 


y 4 = y 3 + t ^(k 1 +k 2 ) = 1.8159 + ^-(1.8159 + 2.1791) = 2.2154 

Càlculo de y 5 

Para n = 4 

k = f{x 4 ,y 4 ) -» k 1 =y 4 =2.2154 

k 2 = f(x 4 + h, y 4 +hk 1 ) = (y 4 + hk^ = 2.2154 + 0.2(2.2154) 

= 2.6585 

y 5 =y 4 + l l(k 1 +k 2 ) = 2.2154 + 2.2154 + 2.6585) 

= 2.7028 


Portanto, temos que y(1.0) = y 1 = 2.7028. 

Exibimos na Tabela 6.2 os valores exatos e os valores aproximados obtidos. 


Tabela 6.2 


N 

*i 

Sol. Exata 

Sol. Aprox. 

Erro 



y(*n) 

y n 

\y(* n )-y n \ 

0 

0 

1.0000 

1.0000 

0.0000 

1 

0.2 

1.2214 

1.2200 

0.0014 

2 

0.4 

1.4918 

1.4884 

0.0034 

3 

0.6 

1.8221 

1.8159 

0.0062 

4 

0.8 

2.2255 

2.2154 

0.0101 

5 

1.0 

2.7183 

2.7028 

0.0155 


Como podemos observar na Tabela 6.2, o Mètodo de Euler Aperfeigoado 
apresentou melhores resultados quando comparado com o Mètodo de Euler, isto 
è, os erros sào menores. 

Para obtermos melhores resultados devemos diminuir a amplitude do passo 
h, por exemplo, tome h = 0.01 e resolva o exercicio da mesma forma, corno an¬ 
teriormente, e compare os resultados obtidos. 
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Observagào 

Os Métodos de Euler e Euler Aperfeigoado sào chamados de Métodos de 
Runge-Kuttade 1 a ordem e 2 a ordem, respectivamente. Esses métodos sào bas¬ 
tante utilizados, pela sua simplicidade e pela precisào dos resultados. 

Existem também, os Métodos de Runge-Kutta de 3 a e 4 a ordem, os quais 
podem ser encontrados com detalhes na Referència [2], 


6.5 Solugào Numèrica de Equagoes Diferenciais Parciais: 
Urna Introdugào 


Como no caso das equagoes diferenciais ordinàrias, vàrios problemas nas 
diversas àreas cientificas podem ser formulados por urna equagào diferencial 
parcial. 


Em geral, as equagoes diferenciais parciais sào mais complexas do que as 
equagoes diferenciais ordinàrias, pois envolvem um nùmero maior de variàveis. 

Assim, quando a fungào incògnita depende de mais de urna variàvel e 
relaciona com suas derivadas parciais através de urna equagào, temos urna 
equagào diferencial a derivadas parciais. 

Por exemplo, a equagào — -2^-^ = 0 è urna equagào diferencial parcial 

dt dx 

do tipo parabòlica em que a fungào incògnita é u = u(x,t). 

Nesta segào, veremos um mètodo numèrico para resolugào das equagoes 
diferenciais parciais, usando para isto urna discretizagào por diferengas finitas. 


Apresentaremos um mètodo de diferengas finitas para resolugào de urna 
equagào diferencial parcial, o qual é baseado na discretizagào do dominio e a 
substituigào das derivadas presentes na equagào diferencial por aproximagòes 
envolvendo valores numéricos da fungào. 


Para isso, inicialmente tratamos a aproximagào das derivadas de urna fungào 
de urna variàvel por diferengas finitas, usando a expansào da Sèrie de Taylor. 

Considere urna fungào y(x) continua com derivadas continuas, podemos 
escrever a fungào y(x) nas vizinhangas de um ponto xpor: 


y(x + h) = y(x) + hy '(x) + 


fi* 

2 


y"(x) + ... 


Ou ainda, 


h 2 


y(x - h) = y(x) - hy\x) + —y"(x) + ... 
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Tornando a expansào da fungào y(x) por Sèrie de Taylor, nas vizinhangas 

do ponto xaté a primeira ordem, e desprezando os termos de ordem 0(h 2 ) e 
de ordem maior, temos: 


y(x + h) = y(x) + hy '(x) 


Dessa forma temos a aproximagào da derivada da fungào y(x), da seguirne 
forma: 


y'(x) 


y(x + h)-y(x) 
h 


Considerando Ay(x) = y(x + h) - y(x), temos a diferenga finita progressiva 

de y(x) e a Fòrmula Progressiva por: 

y(x) = ~ A y(x) 

h 

Analogamente, podemos obter a Fòrmula Regressiva do seguirne modo: 

y'(x)= yix) - yix - h) 
y h 

Considerando Vy(x) = y(x) - y(x - h ), obtemos a diferenga finita regressi¬ 
va de y(x) e a Fòrmula Regressiva por: 

y(x) - T-Vy(x) 
h 

Podemos também, usando a expansào da fungào y(x) por Sèrie de Taylor 
nas vizinhangas do ponto x, até ordem 1, desprezando os termos da ordem 
0(h 2 ) e subtraindo y(x - h) de y(x + h), obter: 

, (x) y(x + h)-y(x-h) 
y 2h 

Considerando 5y(x) = y(x + h) - y(x - h) , temos a diferenga finita centrai 
de y(x) e a Fòrmula Central dada por: 







/'(*) = 


8y(x) 

2h 


Como anteriormente, podemos obter urna aproximagào para y "(x), tornan¬ 
do a expansào da fungào y(x) por Sèrie de Taylor até a ordem 3 e considerando 
a soma das expressòes para y(x + h) e y(x - h ), temos: 


y"(x) 


y(x + h)~ 2y(x) + y(x - h) 
h 2 




Observagào 

As notagòes Ay(x), Vy(x) e 5y(x) representam os operadores de diferenga 
finita progressivo, regressivo e centrai, respectivamente, aplicados à fungào y(x). 

Essas fórmulas de diferengas finitas desenvolvidas para aproximagào das 
derivadas de urna fungào de urna variàvel podem ser utilizadas para aproxima¬ 
gào das derivadas parciais de urna fungào de vàrias variàveis. 

Seja u(x,t) urna fungào de duas variàveis (x,t), temos as seguintes 
diferengas: 

Fòrmula Progressiva na variàvel t 


A A 


Fòrmula Regressiva na variàvel t 


„„u(x,t)-u(x,t-k)_ 

U t * / — / — 



U(x, t ) 


Fòrmula Central na variàvel x 


u(x + h,t)- u(x - h, t) = J_ 

* 2h 2h x 


U xx( X ’t) = 


u(x + h,t)-2u(x,t) + u(x - h,t) 1 


tf 


= -^5 x u(x, t) 










6.5.1 Mètodo de Diferengas Finitas 


Seja urna fungào de duas variàveis u(x,t), considere o plano subdividido 
em retàngulos iguais de lados A x = h e A t = k. 

Considere um ponto P = (x, t) na malha, entào este sera representado por 
P = (x, t) = (x„ tj ) = (ih, jk ). Os indices i e j sào inteiros e os valores h e k sào os 
espagamentos da malha nas diregòes dos eixos x e t. 

Nesse caso, o valor da fungào u(x,t) no ponto P sera denotado por 
u P = u(x j , tj) = u(ih,jk) = u,j. 

O Mètodo das Diferengas Finitas consiste em definir urna malha na regiào 
de interesse a ser calculada, os valores da equagào diferencial parcial e utilizar 
aproximagòes das derivadas vistas anteriormente nessa equagào, num ponto 
genèrico P = (x j ,t j ), conforme mostrado na Figura 6.5. 



Figura 6.5 Malha dos pontos no plano. 

O Mètodo de Diferengas Finitas para urna equagào do tipo parabòlica è 
obtido pelo desenvolvimento da Sèrie de Taylor da fungào u(x, t ) na variàvel t 
(tempo) para formar o quociente de diferengas: 


du 




L/(X;,f +/( )-iV(X;,f.) kd 2 U , 

- k - + -2W {X ^ 


( 2 ) 


para algum e (tj,t j+1 ) e a Sèrie de Taylor para fungào u(x,t) na variàvel 
espacial x: 
























( 3 ) 



emque ^ e(x,_„ x i+1 ). 

Urna vez substituidas as expressòes das derivadas (2) e (3) e desprezando 
os termos de ordem O (k) e de 0(h 2 ) na equagào diferencial parcial, teremos 
urna equagào de diferengas para ser resolvida, escrevendo as derivadas par- 


ciais no ponto P = (x.,f y ) corno: 



II _ Oli _L II 



6.5.2 Mètodo Explicito 


Descreveremos a seguir um mètodo explicito para resolugào de urna equa¬ 
gào diferencial parcial, baseado em diferengas finitas. 

As equagòes diferenciais parciais de segunda ordem classificam-se em 
hiperbólicas, elipticas e parabólicas. Para facilidade de entendimento exibiremos 
um problema unidimensional de difusào dado pela equagào diferencial parabòli¬ 
ca para urna fungào u(x, t ) definida no retàngulo [0, L] x [0, T]: 



(4) 


com condigòes de contorno u(0, t) = u{L, t) = 0 e 0 < t < T, as quais representam 
os valores da fungào u nas extremidades x = 0 e x = L, as condigòes iniciais no 
instante t = 0, u(x, 0) = f(x), 0<x<L, e a é o coeficiente de difusào. 

O Mètodo das Diferengas Finitas Explicito consiste em determinarmos urna 
aproximagào da solugào da equagào diferencial parcial envolvendo urna fungào 
u(x, t ) utilizando, para isto, a derivada parcial de u com relagào a variàvel t (tempo) 
pela diferenga progressiva e a derivada parcial de u com relagào a variàvel x (espa- 
go) pela diferenga centrai, conforme malha da Figura 6.6: 
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Figura 6.6 Malha do Mètodo Explicito. 


Inicialmente, construfmos urna malha no dominio selecionando um nùmero 
inteiro m>0 e definimos um passo nadiregào da variàvel x de tamanho h = l/m. 
Em seguida, selecionamos um passo na diregào da variàvel t de tamanho k. 

Os pontos da malha nesse caso sào do tipo (x ; , f y ), em que x, = ih para 


i = 0,1,...,m, etj=jk, para j = 0,1,... 


A equagào parabòlica (4) implica que para todos os pontos interiores à ma¬ 
lha (. Xj,tj ), paratodo i = 1,2,..., m-1 e j = 1,2,..., teremos: 


du 

~dt 


(Xi,tj) = a 2 



CM;) 


Entào, substituindo as expressòes (2) e (3) em (4) temos a fungào u{x,t) 
nos pontos (x i5 tj), ou seja,t/ (y =t/(x,.,f.) satisfazendo o sistema de equagòes 
lineares corno segue: 


Uj,j + 1 - U iJ =a 2 U J^i- 2U i,j + U i-1,j 


h £ 


(5) 


Para todo i = 1, 2,..., m-1 e j = 1,2,... 

Essa aproximagào leva a um erro de truncamento locai no ponto (/, j) para 
està equagào de diferengas que é dado por: 




















A equagào de difusào discretizada corno em (5) pode ser resolvida para 


u ij+1 conhecendo-se os valores u Lj , u j+1j e u j _ 1j , corno segue: 



( 6 ) 


Para todo / = 1, 2,..., m -1 e j = 1,2,... 




(7) 


Se conhecermos urna condicào inicial u{x, 0 ) = f{x ), para todo 0 < x < 1, te- 
mosque u i0 =f(x i ), para todo i = 0,1,..., m-1 e as condigòes de contorno nas 
extremidades u(0, t) = 0 e u{1, t) = 0 temos que u 0 j = u m J , respectivamente. 

Dessa forma, podemos determinar todos os elementos da forma u h1 ex- 
plicitamente conhecendo as aproximagòes u j0 e continuando o procedimento 
obtemos de maneira similar os valores u j2 ,u i3 ,... 

Exemplo 6.5 

Considere o problema de difusào com valor inicial e condigòes de contorno 
dado pela equagào unidimensional: 



( 8 ) 


dt dx 


Com as condigòes de contorno u(0,t) = u(1,t) = 0,t>0, e com as condi- 
gòes iniciais u(x,0) = sen(nx), 0<x<1. 

Resolvendo numericamente e utilizando o mètodo de diferengas progres- 
sivas na variàvel t e centradas em x, considerando os passos na variàvel x com 
h = 0.2 e com passo na variàvel t com k = 0.01. Observe que a relagào entre os 
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Considerando esses passos temos: 


x j = 0.2i, i = 0,1,2,..., 5 

tj = 0.01j, j = 0,1,2,... 


As condigòes de contorno do problema nos fornecem u 0] = u Xj = 0 e os da- 
dos iniciais, u i0 = sen(n(ih)). Representamos os valores das condigòes iniciais 
naTabela 6.3: 


Tabela 6.3 


i 

j 

x, = ih = 0.2 i 

U U 

0 

0 

0 

0 

1 

0 

0.2 

0.5877852524 

2 

0 

0.4 

0.9510565165 

3 

0 

0,6 

0.9510565163 

4 

0 

0.8 

0.5877852522 

5 

0 

1 

0 


Aplicando a expressào da fòrmula (7): 


u u+ i = 0.25u i+1 j + 0.5u u + 0.25 

ou 

Ui, i+1 = 0-25(u i+1j + 2u i:j + u^j) 

para calcular os valores da solupào aproximada nos pontos da (y' + 7)-ésima 
linha no tempo, conhecendo-se apenas os valores da j-ésima linha. 

Assim, obtemos os elementos u jj+1 para j = 1,2,... e os demais elementos 
a partir das condigòes iniciais e condigòes de contorno. 

Escolhendo inicialmente j = 0, obtemos os valores aproximados da solu- 
gào u 1p u 21 , u 31 e u 41 na primeira linha da malha na variàvel t, corno segue: 

^1,1 ~ 0.25(u 2 o + 2u 10 + u Q 0 ) 

= 0.25(0.9510565165) + 2(0.58778525244 ) + (0) 


= 0.5316567555 














u 21 = 0.25 (u 30 + 2u 20 + u 10 ) = 0.8602387002 


u 31 = 0.25(u 40 + 2u 30 +u 20 ) = 0.8602387005 
u 41 = 0.25{u 5 0 + 2u 4 0 + u 30 ) =0.5316567555 


Da mesma forma, os valores numéricos aproximados de u(x, t ) na segun- 
da linha da variàvel t e x, e em todos os pontos do dominio desejado podem ser 
calculados. 


Como a solucào exata desse problema é dada pela expressào 

u(x,t ) = e ^ f sen(nx), podemos comparar os resultados obtidos numericamen¬ 
te com os valores exatos nos pontos (/, j) e obter os erros correspondentes. 


A estabilidade de um mètodo numèrico està relacionada ao crescimento 
ou decaimento dos erros decorrentes das operagòes aritméticas envolvidas na 
resolugào do sistema de equagòes envolvido. 


Dizemos que um mètodo è estàvel quando esses erros nào aumentam 
quando o nùmero de pontos da malha cresce. 

k 1 

Mostra-se que devemos escolher os valores de h e k de modo que — r < —, 

^ h 2 2 

e jà que a = l no exemplo dado, escolhemos h = 0,2 e k = 0,01 e obtemos 
1 ^ 

\ = ^r = 0.25 <0.5, de modo que essa condigào de estabilidade è satisfeita 
h 

para o mètodo utilizado. 


Caso o valor da relagào entre o espagamento da malha X fosse maior do 
que 0.5, teriamos problemas de estabilidade e buscariamos outro mètodo para 
resolugào. 


6.6 Exercfcios 


1. Usando o Mètodo de Euler, resolva as seguintes equagòes diferenciais 
com valor inicial: 

a) y' = 4-2x, com y(0) = 2. Usando h = 0.01, calcule y(0.05 ). 

b) y' = 1- y/x, com y(2) = 2. Usando h = 0.1, calcule y (2.4). 

c ) y' = -y/x, com y(0 ) = 20. Usando h = 0.2, calcule y(1.0). 




2. Usando o Mètodo de Euler Aperfeigoado, calcule a solugào aproximada 
para as seguintes equagòes diferenciais com valor inicial: 


a) y' = -x/y, com y(0) = 20. Usando h = 0.2, calcule y (1-2). 

b) y' = y-x, com y(0) = 2. Usando h = 0.1, calcule y(0.9). 


3. Com o Software Numèrico - Referència [1] -, usando o Mètodo de Euler, 
calcule a solugào aproximada para as seguintes equagòes diferenciais com valor 
inicial: 

a) y' = -2y + 1, com y(0) = 1. 

Usando h = 0.1 e depois h = 0.01, calcule y(0.9). O que vocè pode afir¬ 
mar sobre a qualidade dos resultados obtidos? Justifique suas afirmagòes 
teoricamente. 

b) y' = y, com y(0) = 1. 

Usando h = 0.2, calcule y(2.0) e um limitante superior para o erro. 

4. Escreva as equagòes discretizadas da equagào diferencial parabòlica da 
difusào para: 

a ) h = 0.1 e k = 0.1 

b)h = 0.5 e k = 0.3 

5. Discretize a equagào diferencial parcial parabòlica da difusào, 


du 

~dt 


(x,t) 


d 2 U 

dx 2 


(X,t) 


utilizando o Mètodo das Diferengas Finitas para valores de passos no tempo 
k = 0.05 e no espago h = 0.1. 

Considere as condigòes de contorno u(0,t) = u(2,t) = 0, 0 <t <T e as con- 
digòes iniciais u(x,0) = 1,0, 0<x<2 para calcular o valor de u(x,t) nos pontos 
(0.9, 0.05), (1,0.05), (1.1,0.05) e (1,0.1). 


6. Discretize a equagào diferencial parcial parabòlica da difusào, 


du 

Ut 


(x,t) = a 2 


d 2 u 

dx 2 


(X,t) 
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utilizando o Mètodo das Diferengas Finitas para valores do coeficiente de difusào 
a-1. 






Escolha valores para os passos no tempo k e no espago h e encontre um 
valor aproximado para a solugào numa regiào quadriculada de lado 1, conside¬ 
rando as condigòes de contorno u(0,t) = u(1,t) = 0, 0 <t <T e as condigòes ini- 
ciais u(x,0) = 1,<x<1. Piote os valores das aproximagòes obtidas nos pontos 
da malha. 

7. Faga um Mapa Conceitual detalhado sobre métodos numéricos para as 
Equagòes Diferenciais Ordinàrias, introduzindo labels indicando se aprendeu 
(A), nào aprendeu (N), se gostou (G), detestou (D ), achou interessante (I), etc. 

8. Faga um Mapa Conceitual detalhado sobre métodos numéricos para as 
Equagòes Diferenciais Parciais, introduzindo labels indicando se aprendeu (A), 
nào aprendeu (N), se gostou (G), detestou (D ), achou interessante (I), etc. 
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